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ABSTRACT 

We present the first results of our deep survey of lensing clusters aimed at constraining the abundance of star-forming galaxies 
at 2 ~ 6 — 10, using lensing magnification to improve the search efficiency and subsequent spectroscopic studies. Deep near-IR 
photometry of two lensing clusters (A1835 and AC114) was obtained with ISAAC/VLT. These images, combined with existing 
data in the optical bands including HST images, were used to select very high redshift candidates at 2 > 6 among the optical- 
dropouts. Photometric selection criteria have been defined based on the well-proven dropout technique, specifically tuned to 
target star-forming galaxies in this redshift domain. 

We have identified 18(8) first and second-category optical dropouts in A1835 (AC114), detected in more than one filter up to 
H (Vega) ~ 23.8 (AB ~ 25.2, uncorrected for lensing). Among them, 8(5) exhibit homogeneous SEDs compatible with star- 
forming galaxies at 2; > 6, and 5(1) are more likely intermediate-redshift EROs based on luminosity considerations. We have also 
identified a number of fainter sources in these fields fulfilling our photometric selection and located around the critical lines. We 
use all these data to make a first attempt at constraining the density of star-forming galaxies present at 6 < z < 10 using lensing 
clusters. Magnification effects and sample incompleteness are addressed through a careful modeling of the lensing clusters. A 
correction was also introduced to account for the expected fraction of false-positive detections among this photometric sample. 
It appears that the number of candidates found in these lensing fields, corrected for magnification, incompleteness and false- 
positive detections, is higher than the one achieved in blank fields with similar photometric depth in the near-IR. The luminosity 
function derived for z > 6 candidates appears compatible with that of LBGs at 2 ~ 3, without any renormalization. The turnover 
observed by Bouwens et al. (2005) towards the bright end relative to the 2 ~ 3 LF is not observed in this sample. Also the 
upper limit for the UV SFR density at 2 ~ 6 — 10, integrated down to I/1500 = 0.3 I/^^s, of p* = 7.4 10~^ Mq yr~^ Mpc~^ is 
compatible with the usual values derived at 2 ~ 5 — 6, but higher than the estimates obtained in the NICMOS Ultra Deep Field 
(UDF). The same holds for the upper limit of the SFR density in the 2 ~ 8 — 10 interval (p* = 1.1 10~^). This systematic trend 
towards the bright end of the LF with respect to blank fields could be due to field-to-field variance, a positive magnification 
bias from intermediate-redshift EROs, and/or residual contamination. Given the low S/N ratio of the high-2 candidates, and 
the large correction factors applied to this sample, increasing the number of blank and lensing fields with ultra-deep near-IR 
photometry is essential to obtain more accurate constraints on the abundance of 2 > 6 galaxies. 
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1. Introduction 

During the last decade considerable advances have been 
made in the exploration of the early Universe, from the 
discovery and detailed studies of redshift z ~ 3 galaxies 
(the so-called Lyman break galaxies, LBGs, e.g. Steidel et 
al. 2003), to z ~ 4-5 galaxies found from numerous deep 
multi-wavelength surveys, to galaxies at z 6-7, close to 
what is believed to be the end of reionization epoch of the 
Universe (e.g. Hu et al. 2002, Kodaira et al. 2003, Cuby et 
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al. 2003, Kneib et al. 2004, Stanway et al. 2004, Bouwens 
et al. 2004b). Extending the searches beyond z ~ 6.5 and 
back to ages where the Universe re-ionized (cf. Fan et al. 
2002) requires extremely deep observations in the near-IR 
bands. Astounding depths can be reached in ultra-deep 
fields, as demonstrated e.g. with J and H imaging of the 
NICMOS Ultra-Deep Field (UDF; Thompson et al. 2005; 
Bouwens et al. 2004a) from which 5 faint {Hab ^ 27) 
candidates at z ~ 7-8 have been identified (Bouwens et 
al. 2004b). 

On the other hand, recent WMAP results seem to place 
the first building blocks of the Universe at redshifts up to 
z ~ 10 — 15 (Spergel et al. 2006) and very distant star- 
forming systems could have been responsible for a signifi- 
cant part of the cosmic reionization. The end of this epoch 
is suggested to be at z ~ 6.0 — 6.5 from the spectrum of 
high redshift quasars (Fan et al. 2002). Therefore, con- 
straining the abundance of z > 7 sources is an important 
challenge of modern cosmology. 

Photometric selection of high-redshift galaxies, such as 
the Lyman-break technique (e.g. Steidel et al. 1995), has 
been successful in identifying star-forming objects at z 
2-4 (cf. Steidel et al. 2003, Shapley et al. 2003) and up to 
z ~ 6 (Bunker et al. 2003, Mobasher et al. 2005). At red- 
shifts z > 5, only w 30 galaxies are currently known with 
confirmed redshifts (cf. review by Spinrad 2003). For now, 
more than 5 galaxies with z ^ 6.5 (Hu et al. 2002, Kodaira 
et al. 2003, Cuby et al. 2003, Kneib et al. 2004) are known 
with secure redshifts. These objects are generally found 
through their Lyman-a emission which is redshifted into 
the ~ 9200 A window, the "reddest" window relatively 
free of skylines in the optical. The abundance of z '--^ 10 
galaxies was recently discussed by Bouwens et al. (2005a) 
using NICMOS-UDF Jno and TJieo data. Their conclu- 
sion is that strong evolution exists between z ~ 7 — 8 and 
z ~ 3 — 4, the SFR density being much lower at very high 
z down to the limits of their survey (iisoo — 0.3L*^3). 
However, it is crucial to increase both the size and the 
depth of the surveyed field to set strong constraints on 
the star-formation at z > 7, as field-to-field variance can 
be important. 

In this paper we present the first results of our deep 
survey of lensing clusters aimed at constraining the abun- 
dance of star- forming galaxies at z 6 — 10, using lensing 
magnification to improve the search efficiency and subse- 
quent spectroscopic studies. The motivations are the fol- 
lowing. On the one hand, our understanding of the first 
generation of stars and galaxies, the so-called Population 
III objects (cf. review of Loeb & Barkana 2001), has im- 
proved with the development of new models for these low- 
metallicity starbursts (Tumlinson & ShuU 2000, Bromm et 
al. 2001, Schaerer 2002, Schaerer 2003). The observational 
properties inferred from this modeling show us that it is 
now possible to detect some of these objects at z ~ 8 — 10, 
thanks to the intensity of their emission lines, using very 
deep near-IR surveys on 8 m-class telescopes (e.g. Schaerer 
& Pello 2001, Barton et al. 2004). 



Our project is based on the photometric pre-selection 
of candidates making use of the natural magnification pro- 
vided by foreground lensing clusters. This technique, first 
referred to as the "gravitational telescope" by Zwicky, has 
proven highly successful in identifying a large fraction of 
the most distant galaxies known today thanks to magni- 
fications by typically 1-3 magnitudes (e.g. Franx et al. 
1997, Elhs et al. 2001, Hu et al. 2002, Kneib et al. 2004). 

We present a color-color selection similar to the 
Lyman-Break technique, used to identify very high red- 
shift objects using their specific spectrophotometric prop- 
erties. As first targets for our survey, we have chosen fields 
centered on lensing clusters with well-constrained mass 
distributions, and already known to be efficient gravita- 
tional telescopes. We use the sample of high-redshift can- 
didates selected in these fields to constrain the abundance 
of star-forming galaxies up to z ~ 10. 

The plan of this paper is as follows: in Sect. El we 
justify the observing strategy and the photometric crite- 
ria adopted in this project. Photometric observations and 
data reduction are described in detail in Sect. O and 2| 
The construction and analysis of the photometric catalogs 
is given in Sect.El The photometric selection of very high 
redshift candidates is presented in Sect. El The properties 
of the final list of candidates, including spectral energy 
distributions (hereafter SEDs), photometric redshifts and 
magnification estimates, are detailed in Sect. [3 In Sect. 
(HI we discuss the intrinsic (lens-corrected) properties of 
this sample, the number-density of star-forming galaxies 
at z ~ 6 — 10 as compared to simple model expectations, 
and the implications for the cosmic SFR. We also sum- 
marize our ongoing spectroscopic survey in these fields. 
Conclusions are given in Sect. El In the appendix wc pro- 
vide more details on the improvement of the data reduc- 
tion procedure, the completeness and false-positive detec- 
tion estimates in the different fields and filters, as well 
as additional tests performed on the reliability of optical 
dropouts. 

Throughout this paper, we adopt the following cosmol- 
ogy: a flat A-dominated Universe with the values Q\ = 
0.7, n„, = 0.3, Qb = 0.045, Hq = 70 km Mpc'^, and 
(Tg = 0.9. All magnitudes given in the paper are quoted 
in the Vega system. Conversion values between Vega and 
AB systems for these filters are given in Table 

2. Simulations and observing strategy 

Our project aims to search for z > 6 galaxies, typically 
in the z ~ 7-10 domain. We have performed simula- 
tions to estimate the expected magnitudes of galaxies at 
such redshifts, and to establish robust photometric criteria 
to select high-redshift candidates behind lensing clusters. 
For this, we have used the evolutionary synthesis star- 
burst models by Schaerer (2002, 2003) for Population III 
and extremely metal deficient galaxies, together with the 
usual templates available for normal galaxies. In particu- 
lar, we used the empirical SEDs compiled by Coleman, 
Wu and Weedman (1980) to represent the local popu- 
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lation of galaxies, with spectra extended to wavelengths 
A < 1400 A and A > 10000 A using the equivalent spec- 
tra from the Bruzual & Chariot GISSEL library for solar 
metallicity (Bruzual & Chariot 1993). We also included 
the starbursts templates SBl and SB2, from Kinney et 
al. (1996), and the low metallicity galaxy SBS0335-052 
(Izotov, 2001, private communication). 

We consider a fiducial stellar mass halo of 10^ Mq, cor- 
responding to a collapsing DM halo of 2 x 10^ Mq. With 
a ACDM model and the cosmological parameters adopted 
here, we expect ^ 10 of such DM haloes per Mpc'^ within 
the relevant redshift range z ~ 5 — 10, and ~ 1 Mpc~'^ 
with a DM halo of 10^ Mq (e.g. Loeb & Barkana 2001). 
The virial radius is of the order of a few kpc, and thus we 
consider that sources are unresolved on a 0.3" scale, with 
spherical symmetry. The reionization redshift is assumed 
to be z ~ 6, but results discussed below are independent 
of this precise value. Lyman series troughs (Hainian & 
Loeb 1999), and Lyman forest absorptions following the 
prescription of Madau (1995) are included. Although a de- 
tailed description of Lyman-a emission is out of the scope 
of this paper, we have taken into account the possible im- 
pact of the main emission lines on the integrated colors us- 
ing rather simple and extreme assumptions. Simulations 
accounting for an extended Lyman-a halo (cf. Loeb & 
Rybicki 1999) have been computed, together with a simple 
parametrization of the fraction of Lyman-a flux entering 
the integration aperture. Two extreme assumptions are 
considered here for the IMF, either a standard Salpeter 
IMF, with stars forming between 1 and 100 Mq, or a top- 
heavy IMF, with stellar masses ranging between 50 and 
500 Mq. Some preliminary versions of these simulations 
were presented in Schaerer & Pello (2001) and Pello & 
Schaerer (2002). Figure ^ shows, for each IMF, the ex- 
pected SED of a 10^ Mq stellar mass halo, using two dif- 
ferent assumptions for age. 

Nebular continuous emission and strong emission lines 
could have important effects on the integrated fluxes and 
colors of such galaxies, although broad-band colors alone 
do not allow one to precisely constrain the physical prop- 
erties of these galaxies. The main signatures of genuine 
star-forming sources at z > 6 are common to all models: 
they are optical dropouts, displaying a strong break and 
"red" optical vs. IR colors, whereas they exhibit a "blue" 
SED longward of Ly-a provided reddening is sufficiently 
small. Different redshift intervals can be defined using an 
appropriate set of near-IR filters in combination with opti- 
cal data. This particular application of the Lyman break 
technique as a function of redshift is shown in Figs. |21 
to 01 for the different redshift intervals considered in this 
paper. For clarity, only a reduced number of models is 
presented in these figures. Color shifts corresponding to 
an intrinsic extinction of Ay — 1 (Calzetti et al. 2000) 
are shown by arrows, for starbursts at z = 1 and z = 3. 
Stellar colors presented in these diagrams were computed 
from the library of Pickles (1998). 

Figure |21 displays the J — H versus H — Ks color-color 
diagram for different extreme Population III starbursts 
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Fig. 1. Examples of SED for Pop III models used in Fig. 
Eland Fig. El (upper panel), for a fiducial 10^ Mq burst 
at z = 9, with ages 10* and 10^ yrs (solid and dashed 
lines), and two IMF (normal Salpeter -black- and top- 
heavy -red-), compared to the transmission curves of the 
FORS/ISAAC filters set used in this survey (lower panel). 
The top axis gives the corresponding redshift at the wave- 
length of the Lyman-a break. We overplot the location 
of the Hell 1640 emission line (without scaling its fiux), 
which has a small impact on the k-correction as shown in 
Figs. O and Fig. Different assumptions were considered 
for the Lyman-a emission. 



within the 5 < z < 11 interval, compared to the loca- 
tion of stars and normal galaxies at different redshifts. 
This diagram is particularly well suited to identify galaxy 
candidates in the 8 < z < 11 interval among the opti- 
cal dropouts. At redshifts above ^ 10, galaxies are not 
detected in the J band (see also Fig. [SJ • For galaxies at 
6 < z < 9, the same photometric selection can be per- 
formed including the z (0.9 /im) and SZ (1.1 /im) filters 
(Fig.|21and0J. The characteristics of these filters are sum- 
marized in Tableland Fig. ^ 

Once the color-selection is well established, an impor- 
tant issue is the photometric depth needed to detect typ- 
ical stellar haloes up to a given mass, in order to derive 
statistically significant results on the nature and proper- 
ties of z > 6 sources. 

According to our simulations, the predicted magni- 
tudes in the Vega system for a reference stellar halo mass 
of 10^ Mq, zero metallicity stars, a top-heavy Salpeter 
IMF, and a starburst younger than 10^ yrs, typically range 
between ~ 24.5 and 26.0 in J in the relevant redshift range 
(z < 8), ~ 24.5 and 25.5 in H, and - 24 to 25 in Ks 
(z < 10), depending on models, within the z ~ 6 — 10 in- 
terval (see Fig. Ol. For the same models, SZ and z range 
between ~ 25 — 26, for z < 7 and z < 6 respectively, and 
sources become dropouts in these filters beyond these red- 
shifts (see Fig. Eland Fig.P). For a standard Salpeter IMF 
(from 1-100 A^q) these values are ~ 2 magnitudes fainter 
than for a top-heavy IMF; increasing the metallicity for 
this IMF implies a somewhat larger UV restframe flux (up 
to 0.5 mag brighter for solar metallicity; cf. Schaerer 2003, 
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Fig. 2. J dropout selection in the z ~ 8-11 domain. J — H 
versus H — Ks color-color diagram (Vega system) showing 
the position expected for different objects over the interval 
z to 11. The position of stars and normal galaxies up 
to z < 8 are shown, as well as the shift direction induced 
by Ay — 1 magnitude extinction. Thin and thick lines 
display models below and above z = 8. Several models for 
Pop III starbursts are presented, for different fractions of 
Lyman-a emission flux entering the integration aperture: 
100% (red solid line), 50% (red dashed hue) and 0% (red 
dot-dashed line). The location of Kinney et al. (1996) star- 
bursts templates is also given for comparison (SBl(cyan) 
and SB2 (blue)). All star- forming models enter the high-z 
candidate region at z > 8. 
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Fig. 3. SZ dropout selection in the z ^ 7-8.5 domain. 
SZ — J versus J — H color-color diagram (Vega system) 
showing the position expected for different objects over 
the interval of z ^ to 8.5. Thin and thick lines display 
models below and above z = 7 respectively. Models dis- 
played and general comments are the same as in Fig. |21 
The position of stars and normal galaxies up to z < 7 are 
shown. 



Fig. 2). Also, obviously these magnitudes scale with stel- 
lar mass. In other words, a stellar halo with a standard 
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Fig. 4. z dropout selection in the z ~ 6-7.5 domain, 
z — SZ versus SZ — J color-color diagram (Vega system) 
showing the position expected for different objects over 
the interval of z ^ to 7.5. Thin and thick lines display 
models below and above z = 6 respectively. Models dis- 
played and general comments are the same as in Fig. |5| 
The position of stars and normal galaxies up to z < 6 are 
shown. 



IMF exhibits about the same magnitudes as a top-heavy 
IMF which is a factor of 10 less massive, all the other pa- 
rameters being the same. This comment not only stands 
for PopIII models, but also for solar metallicity starbursts 
and constant star-forming models with standard IMF and 
metallicity. 

Assuming a minimum gravitational magnification of 
~1 magnitude, if we intend to detect stellar haloes up to 
10^ Mq (or a few lO"^ Mq, depending on IMF), the pho- 
tometric depth required is of the order oi H ^ 24.0 and 
Ks ~ 23.5 for a positive detection, and up to ~ 26.0 in z, 
SZ and J to identify significant dropout sources in these 
filters. The number of sources expected can be roughly es- 
timated as follows. Taking into account the typical covol- 
ume surveyed in a lensing cluster under these conditions 
^ , the relevant density of DM haloes assuming a conserva- 
tive fraction of 10% of baryonic mass converted into stars 
before z = 6 (^ 0.1 Mpc^, corresponding to 10^" Mq DM 
halo), and the probability of detection related to the visi- 
bility of the starburst within the relevant redshift interval 
(10^ yrs restframe, thus a duty-cycle factor ~ 0.1 to a few 
0.01; see Sect.|Sl), the number of sources expected ranges 
between a few tens and a few hundreds. We still expect a 
few positive detections at z = 6 — 10 with a completeness 
level of the order of 10% or even lower. However, strong 
lensing effects have to be carefully taken into account in 
this survey, as explained below. A detailed comparison be- 
tween the number of sources expected and the number of 



^ a few 10* Mpc"^ between z — Q and 10, for a ~ 2 x 2 arcmin^ 
field of view, after correction for a typical magnification factor 
of ~ 2 
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sources actually detected in our lensing fields is provided 
in Sect. El 



3. Photometric data 

Two lensing clusters were selected for tfiis pilot study with 
the VLT: AC114 (a=22:58:48.26 (5=-34:48:08.3 J2000, 
z = 0.312) and Abell 1835 (a=14:01:02.08 ^=+02:52:42.9 
J2000, z = 0.252). AC114 is a well-known "gravita- 
tional telescope" , for which multiwavelength observa- 
tions are available. The corresponding lens model is well- 
constrained by a large number of multiple-images spec- 
troscopically identified at high redshift (Smail et al. 1995, 
Natarajan et al. 1998, Campusano et al. 2001, Lemoine- 
BusseroUe et al. 2003). Its Einstein radius is about 30" for 
z > 6. Abell 1835 is the most X-ray luminous cluster in 
the XBACS sample (X-ray- Brightest Abell-type Clusters 
of galaxies, Ebeling et al. 1998), thus potentially one of the 
most efficient gravitational telescopes. Indeed, strongly 
lensed features were identified in this cluster, based on 
deep ground-based and HST images, and it was used to 
search for bright submm galaxies with SCUBA (Smail et 
al. 1999, Ivison et al. 2000). The mass model used is simi- 
lar to the one developed by Smith et al. (2005), and gives 
an Einstein radius of ~ 40" at high z. 

We observed these clusters with ISAAC and FORS 
in the near-infrared domain 0.9 to 2.2 /im) between 
September 2002 and April 2004, covering as far as possi- 
ble the z, SZ, J, H, and K bands. Transmission curves 
for these filters are presented in figure ^ In addition, op- 
tical images at shorter wavelengths (from U to I band) 
are available in our group from previous surveys, or from 
archival data. Table^summarizes the main characteristics 
of the photometric dataset. 

Note that in the remainder of the paper we shall loosely 
use the term "near-lR" for the SZ, J, H , and K filters, 
whereas "optical" refers to all filters shortwards of 1.0 /im, 
from U to I (or z-band when available). 

3.1. Near-IR photometric observations 

We obtained imaging data with the Infrared Spectrometer 
And Array Camera (ISAAC, Moorwood 1997) located on 
the Nasmyth-B focus of the 8.2m VLT-UTl (Antu tele- 
scope), using the Short- Wavelength channel of the instru- 
ment (Cuby et al. 2000). The field of view of the camera is 
about 2.5 arcmin x 2.5 arcmin with a pixel size of 0.148". 
The data for AC114 were acquired during UT 2002 August 
19-20-22 (period 69). Due to technical problems with the 
instrument in this period, the usual ISAAC detector was 
changed for the Aladdin 1024 x 1024 InSb array. The data 
for Abell 1835 were acquired during UT 2003 January 
14, February 9-11-12-14-15 (period 70, JHK), and UT 
2004 April 20 and May 15 (period 73, SZ) with the usual 
Hawau Rockwell 1024 x 1024 Hg:Cd:Te Array. Differences 
in efficiency between these two detectors have been re- 
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Fig. 5. From top to bottom, J, H and Ks magnitudes as a 
function of redshift for a top-heavy IMF, for a fiducial stel- 
lar halo of 10^ Mq. The values corresponding to a normal 
Salpeter IMF are about 2 magnitudes fainter over all the 
redshift interval. Black and red lines correspond respec- 
tively to burst ages 10^ and 10** yrs. Various models for 
Pop 111 objects are presented, for different fractions of the 
Lyman-a emission entering the integration aperture: % 
(thick dot-dashed line), 50 % (thick dashed line), and 100 
% (thick solid line). Thin dot-dashed lines correspond to 
a self-consistent extended Lyman-a halo emission (Loeb 
& Rybicki 1999), whereas thin dashed lines display the 
same model with 100 % of Lyman-a emission entering the 
integration aperture. 
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Fig. 6. From top to bottom, z and SZ magnitudes as a 
function of redshift, for a fiducial stellar halo of 10^ M0 
and the same models as in Fig. O 



3.2. Optical and intermediate band images 

z band observations of Abell 1835 were obtained dur- 
ing UT 2004 March 26 and April 10 with the FOcal 
Reducer/low dispersion Spectrograph (F0RS2) at VLT. 
This instrument has a 0.252" pixel size and a field of view 
of 7.2 arcmin x 7.2 arcmin. We used dithered individual 
exposures of 120 s. 

To be able to distinguish between "low" (z ^ 6) 
and "high" (z > 6) redshift objects, we have compiled 
the available optical images for the two lensing clusters. 
For AC114, we used the data from Campusano et al. 
(2001) covering U to I filters, including a mosaic of deep 
HST/WFPC2-F702W (R) observations. Images in this 
band were obtained in both high-sky and low-sky modes, 
with different orientations on the sky, in such a way that 
the final composite image in this filter exhibits three re- 
gions with different behavior (noise properties, photomet- 
ric depth) across the ISAAC field of view. Each of them 
is presented separately in Tabled 

Optical data for Abell 1835 include BVRI imag- 
ing obtained with the CFH12k camera at CHFT, and 
HST/WFPC2-F702W (i?) images acquired in low sky 
mode. Because of the field centering chosen for the near- 
infrared data, only half of the ISAAC field of view is cov- 
ered by the HST/WFPC2 image. 

For the two clusters, all optical data fully cover the rel- 
evant region studied around the cluster center. The entire 
ISAAC field is covered in most cases. The overlap fraction 
between the optical images and the near-IR data is indi- 
cated in Tabled where the references and main properties 
of the data set are summarized. Overlap fractions refer to 
the ISAAC field. 



ported, the Hawaii Rockwell detector providing better re- 
sults in terms of photometric accuracy. ^ 

Near infrared imaging is challenging because of the 
dominant and variable sky background. We used dithering 
of short exposures with subintegration x integration times 
of 4 X 45 s in the SZ and J bands, 11 x 12 s in the H 
band, and 6 x 15 s in the Ks band, with a 30" jitter 
box. These values provided a good compromise between 
an optimal photometric depth over a large fraction of the 
field of view, and good sky-subtraction in a crowded field. 
For Abell 1835, the field center was chosen such as to 
prevent contamination by a very bright star located at 
the north of the cluster. Therefore the brightest cluster 
galaxy is not at the center of the field (see Fig. llOfl. 

Calibration data were obtained in the usual way (de- 
tector darks, twilight flats, ...). Standard stars from the 
LCO/Palomar NICMOS hst (Persson et al. 1998) were 
used for photometric calibration. 



^ For more details see Section 1.2.2 of the ISAAC user man- 
ual (Cuby et al. 2002) 



4. Data reduction and calibration 

Near infrared photometry of extremely faint sources re- 
quires a careful data reduction. The general procedure de- 
scribed here was performed for all the ISAAC data {SZ , 
J , H and Ks bands). A number of specific improvements 
are given with more details in App. El For the F0RS2 (z 
band) data, we used a standard flat-field correction and 
combination of the individual frames with bad-pixel rejec- 
tion. 

We reduced our data using IRAF procedures and ac- 
cording to the ISAAC Data Reduction Guide v. 1.5^. The 
different steps are the following: photometric calibration, 
bias subtraction, flat-fielding, sky subtraction, registration 
and combination of the images. The reduction recipe we 
used, as well as some of the improvements, were mostly 
inspired from the reduction of near-IR observations on the 
HDFS field with the same instrument (Labbe et al. 2003). 

Photometric zero-points were derived from 
LCO/Palomar NICMOS standard stars (Persson et 
al. 1998), observed each night using a five-point jitter 
pattern. After subtracting from each image the median 

See http : //www . eso . org/ instruments/ Isaac/ 1 
drg/html/ drg . html 
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Table 1. Main properties of the photometric dataset used in this paper: filter identification, total exposure time, 
average seeing measured on the original images, pixel size, 1 a limiting magnitude inside a 1.5 " diameter aperture, 
filter effective wavelength, AB correction, overlap fraction relative to the ISAAC frames (covering an area of 2.5 arcmin 
X 2.5 arcmin), and references. Exposure time values of three different subsets of the HST-_R702vk image, annotated 
i?7Q2'^ are converted into low sky mode for comparison (see text for details). AB corrections {Cab) correspond to 

rriAB = TTlVega + CaB- 



Filter 


iexp 


seeing 


pix 


depth 




Cab 


overlap 


Reference 




[ksec] 


["] 


["] 


[mag] 


[nm] 


[mag] 


[%] 




AC114 


U 


20.00 


1.3 


0.36 


29.1 


365 


0.693 


43.1 


Barger et al. 1996 


B 


9.00 


1.2 


0.39 


29.0 


443 


-0.064 


100.0 


Couch et al. 2001 


V 


21.60 


1.1 


0.47 


28.5 


547 


0.022 


76.6 


Small et al. 1991 


R702 


> 8.30 


0.13 


0.100 


> 27.7 


700 


0.299 


84.9 


Natarajan et al. 1998 


R702 


> 24.90 


0.13 


0.100 


> 28.4 


700 


0.299 


41.4 


Natarajan et al. 1998 


R702 


40.00 


0.13 


0.100 


28.6 


700 


0.299 


17.0 


Natarajan et al. 1998 


Isi4 


20.70 


0.3 


0.100 


26.8 


801 


0.439 


77.6 


Small et al. 1991 


J 


6.48 


0.52 


0.148 


25.5 


1259 


0.945 


100.0 


This work 


H 


13.86 


0.40 


0.148 


24.7 


1656 


1.412 


100.0 


This work 


Ks 


18.99 


0.34 


0.148 


24.3 


2167 


1.873 


100.0 


This work 


Abell 1835 


V 


3.75 


0.76 


0.206 


28.1 


543 


0.018 


100.0 


Czoske et al. 2002 


R 


5.40 


0.69 


0.206 


27.8 


664 


0.246 


100.0 


Czoske et al. 2002 


R702 


7.50 


0.12 


0.100 


27.7 


700 


0.299 


45.7 


Smith et al. 2005 


I 


4.50 


0.78 


0.206 


26.7 


817 


0.462 


100.0 


Czoske et al. 2002 


z 


6.36 


0.70 


0.252 


26.7 


919 


0.554 


100.0 


This work 


SZ 


21.96 


0.54 


0.148 


26.9 


1063 


0.691 


100.0 


This work 


J 


6.48 


0.65 


0.148 


25.6 


1259 


0.945 


100.0 


This work 


H 


13.86 


0.50 


0.148 


24.7 


1656 


1.412 


100.0 


This work 


Ks 


18.99 


0.38 


0.148 


24.7 


2167 


1.873 


100.0 


This work 



sky of all the pattern, we measured the total counts 
in a 20 pixel radius circular aperture, and from these 
integrated fluxes we derived the zero-points. Airmass 
differences between science exposures were corrected to 
a reference value for each filter using a linear relation 
between zero-point and airmass derived from standards 
stars observed at different airmasses. 

After removing the instrumental ghost, substracting 
a median dark frame and flat-fielding our data, we used 
the IRAF package XDIMSUM** to apply a two-step sky- 
subtraction. During the first pass, each image is sky- 
subtracted using the sky pattern obtained from a group of 
adjacent frames and a bad pixel mask is created in the pro- 
cess. The relative shifts between images are derived from 
the position of several stars matched in each frame. Then, 
images are registered and combined using integer shifts 
values to preserve the noise properties and rejecting all 
bad pixels. Bright sources are detected in order to create 
an object mask, and a second sky-subtraction is applied to 
the data, this time using the mask to reject pixels located 
on objects in the evaluation of the sky. This improves the 
quality of the final stacked image. 

As a cross-check for our stacking procedure, we pro- 
duced another version of the final images using a slightly 



* XDIMSUM is a modified version by the IRAF group of the 
Deep Infrared Mosaicing Software package by P. Eisenhardt et 
al. See ,ftp: / / iraf.noao.edu / extern-v212/ xdimsum, for details 



different reduction recipe. After the usual ghost and dark 
removals, we flat-fielded each image using a sky flat, cre- 
ated by evaluating the sky in a group of adjacent frames, 
again masking the bright objects. Then individual frames 
were registered and coadded in a standard way. The re- 
sulting images are found to be similar, in terms of qual- 
ity, detection level and photometric depth, to the general 
procedure described above. The main difference is an en- 
hanced quality around the bright galaxy haloes close to 
the cluster core, and thus we finally adopted the two-step 
sky-subtraction procedure. 



5. Analysis of images 

Since we use imaging data acquired with very different 
filters and instruments, we had to match them to a com- 
mon reference when measuring the required multi-band 
photometry. However, the registration and seeing match- 
ing process generates the resampling of data, and there- 
fore modifies the noise properties of the background in the 
sense that the error bars measured on these modified im- 
ages by the standard means tend to be underestimated. 
For this reason, we preferred to use the original images to 
derive the error bars in each band, as explained below. 
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5.1. Image registration and astrometry 

We registered the final SZ, J, iJ,and Ks band images 
for each cluster with a simple shift, except in the case 
of AC114, where we corrected for a slight distortion that 
appeared in the J band image. 

In order to measure relative photometry in the near- 
IR bands, we matched all the images to a common seeing 
using a simple Gaussian convolution, the worst case being 
the J band for both clusters. 

As all the photometric high-z candidates are expected 
to be detected in the H band we have defined an H band 
selected sample. The detection images were created with 
the original H band images, weighted by the square root 
of the corresponding exposure time maps in order to get 
a uniform background noise across the field. 

The available optical images for each cluster were reg- 
istered to the ISAAC combined images, using standard 
IRAF procedures for rotation, magnification and resam- 
pling of the data. These images were mainly used to ex- 
clude well-detected low-z sources, and also for the first 
visual inspection of the optical dropouts. However, resam- 
pling of data could produce both spurious detections and 
false non-detections among the faint sources. For this rea- 
son, further discussed in Sect. 16.21 we have used the orig- 
inal images instead to define our final sample of optical 
dropouts and we work on object coordinates. 

We performed an astrometric calibration for all these 
images using ~ 30 bright unsaturated objects present in 
the USNO catalog (USNO A2.0, Monet et al. 1998). The 
error obtained in the absolute astrometric calibration is 
typically 0.2 " for a whole ISAAC field of view. 

5.2. Photometry 

We used the S Extractor package version 2.2.2 (Berlin 
& Arnouts 1996) to detect objects and to compute mag- 
nitudes within our images. We optimized the parame- 
ters to detect very faint unresolved sources. Magnitudes 
were measured within identical circular apertures in all fil- 
ters (from optical to near-infrared), with the S Extractor 
"double-image" mode, using the H band detection im- 
ages described in Sect. 15.11 Near infrared images were 
seeing-matched to the J band (0.52"or 0.65", see Tablenj, 
but the seeing measured on ground-based optical bands is 
worse. Since optical images were only used here for non- 
detection purposes (i.e., for the identification of optical 
dropouts), we preferred to keep their original seeing in- 
stead of degrading the quality of near-IR and HST images. 

S Extractor detection parameters were the following: 4 
connected pixels above a threshold of 1 cr on the detection 
image, which corresponds to a central value of about 4.5 a 
for a seeing-limited source, where a stands for the typical 
local background noise. Magnitudes were measured within 
a 1.5"-diameter aperture (i.e., 10 pixels on the ISAAC im- 
ages). This value is an optimized compromise allowing us 
to obtain a fair estimate of the total magnitude for point- 
like sources on near-IR images, while keeping a good S/N 



on the optical images with the worst seeing (r-^ 1"). We 
have also checked that magnitudes and error-bars mea- 
sured within a 12 pixels (1.7") diameter aperture are usu- 
ally consistent within Icr error bars for the faintest sources 
considered in this study. 

Since SExtractor uses the registered, seeing-matched 
images to compute the photometric errors in all bands, 
these values are systematically too optimistic, thus lead- 
ing to artificially high S/N determinations. To get more 
realistic errors in our photometry, we preferred an empiri- 
cal method to derive them. We measured the typical RMS 
in the pixel distribution within apertures of the same phys- 
ical size as for flux measurements: we averaged the pixel to 
pixel statistics in about 5000 non-overlapping apertures, 
randomly thrown inside regions free of objects in each orig- 
inal (unregistered, unconvolved) image. For each entry in 
the SExtractor output catalog, we recomputed the pho- 
tometric error using the value of a derived from these sim- 
ulations. This photometric error measured in the original 
images was also used to compute the limiting magnitude 
in each band, reported in TableQ] All S/N values reported 
throughout the paper refer to these mock simulations. 

The final catalogs include photometry within 1.5" 
aperture for all objects detected in the H band; we were 
able to measure photometry of very faint sources (J ^ 
24.4 - 24.8, H and Ks - 23.5) with a relatively good 
accuracy (S /N > 3 — 4) . The effective exposure time to- 
wards the edges of the field is smaller than in the cen- 
tral region due to the dithering pattern used in near-IR 
images, thus leading to brighter detection limits and an 
increasing number of spurious detections at the edges of 
the frame. In this study, we use only the region of the field 
for which the effective exposure time is above 50 % of the 
(maximum) total value. The overlap fractions relative to 
the ISAAC frames given in Table prefer to these central 
regions, corresponding to 6.34 arcmin^ for Abell 1835 and 
6.10 arcmin^ for AC 114. The images shown in Fig. 1101 and 
^2 also refer to these central regions. 

5.3. Photometric completeness 

The characteristics of the final processed images are pre- 
sented in Table Q The seeing was measured on the origi- 
nal co-added images. We have computed the completeness 
values for point-sources, in each cluster and band, for near- 
IR magnitudes within the relevant intervals. These limits 
were obtained from direct simulations as follows. Artificial 
stars (i.e., seeing limited sources) of fixed magnitude, rang- 
ing between 21 and 25, were added 1000 times at 30 differ- 
ent random locations on our images, and then extracted 
using the same method for detection and photometry as 
for astronomical sources (described above). Only "free" 
sky regions were used for this exercise, with uniform noise 
properties. This excludes in particular the cluster core and 
the edges of the images, where the effective exposure time 
is less than 50 %. Completeness values are derived from 
the fraction of objects we recovered in our images. The 
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corresponding curves are shown in Fig. 13 and the com- 
pleteness levels are reported in Tabic El A completeness 
level of 20 — 30% is achieved for H and Ks ~ 23, thus 
in good agreement with the requirements given in Sect.|21 
for at least a few positive detections at z ~ 6 — 10. 

Table 2. 90 % and 50 % completeness limits for each 
cluster and band, corresponding to simulations shown in 
Fig.CI 



Cluster 


Band 


90 % compl. 


50 % compl. 


AC114 


J 


23.5 


23.9 


AC114 


H 


22.5 


22.8 


AC114 


Ks 


22.4 


22.8 


Abell 1835 


z 


23.0 


23.6 


Abell 1835 


sz 


24.0 


24.4 


Abell 1835 


J 


23.3 


23.6 


Abell 1835 


H 


22.3 


22.7 


Abell 1835 


Ks 


22.1 


22.7 




Magnitude (Vega system) 

Fig. 7. Completeness curves obtained from simulations in 
each cluster, for the z band (orange diamonds), SZ band 
(blue stars), J band (black triangles), H band (red circles) 
and Ks band (green squares). 

Our near-IR survey has reached SZ ~ 25.6, J ^ 24.3- 
24.4, H - 23.5 and Ks - 23.1 (AC114) to 23.5 (AbcU 
1835) (3 a detection level within 1.5" aperture), i.e. AB~ 
25-25.5 in JHKs and AB~26.3 in SZ. The minimum 
magnification factor over the region covered by our near- 
IR survey is ~ 0.7 magnitudes, and at least ~ 1 magnitude 



over 50 % of the ISAAC field of view. Thus, the effective 
3cr limiting magnitudes reached here are close or similar to 
those attained in the HDFS (Labbe et al. 2003) in JHKs 
(respectively AB~ 26.8, 26.2 and 26.2). Note that the lim- 
iting Vega magnitudes given in Table ^ correspond to la- 
values. Our 3(T limiting magnitudes in the H band also 
comeclose to the magnitudes of the z ~ 7 — 8 z-dropouts 
detected by Bouwens et al. (2004b) in the Hubble Ultra 
Deep Field, with Hiqo{AB) ~ 26.0 to 27.3, after correction 
for a typical magnification factor of at least 1 magni- 
tude. 

We determine below the additional correction for in- 
completeness generated by our sample selection proce- 
dure. Lensing models were used to derive the effective 
completeness of our sample as a function of redshift and 
magnitude, as compared to blank fields. 

6. Selection of high-z photometric candidates 

This Section presents the procedure adopted to select the 
sample of high-redshift candidates. As described in the 
observing strategy (Sect.|21), color-color diagrams of opti- 
cal dropouts have been used to select high-redshift can- 
didates from our deep near-infrared images. For a sub- 
sample of them, individual probability distributions and 
photometric redshifts can be reliably derived from their 
photometric SEDs, as discussed in Sect. [3 The location 
of the photometric candidates with respect to the critical 
lines, and thus the typical magnification factors reached 
by our sample, are also briefly described. 

6.1. Near-IR color-color diagrams of bright objects 

As an additional test of the photometric selection of 
sources, we have checked that magnitudes and colors of 
bright sources {H < 22.5, i.e. S/N > 8 — 10) are in good 
agreement with expectations. We have already secured the 
colors of cluster members during the photometric calibra- 
tion described in Appendix^ Stars morphologically iden- 
tified by SExtractor {flag^, > 0.95) are located at the 
expected position in these diagrams. A representative ex- 
ample is given in Fig.lHlfor AbeU 1835. This diagram is to 
be compared to the theoretical expectations displayed in 

Fig.ia 

If we adopt the {R — Ks > 5.6) definition from Daddi 
et al. (2000) to select Extremely Red Objects (hereafter 
EROs) from this sample of bright sources, we find that 
a few of them lie inside the low-redshift region of the 
diagram for each cluster, a location mainly compatible 
with dust-reddened starbursts. This diagram also shows 
that only very few bright objects correspond to our color- 
selection criteria for high redshift galaxies. All these ob- 
jects are detected in the optical bands. We have carefully 
inspected the morphology of these sources on the HST im- 
ages and found two cases: either they clearly correspond 
to two blended objects, or they are point-like sources. In 
the first case, the two objects merge in a single source on 
the near-IR images, in which case the photometric mea- 
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Fig. 8. Location of objects brighter than H = 22.5 in 
the J — H vs H — Ks color-color diagram for Abell 1835 
(black dots). Special symbols are used for S Extractor 
stars {flag,: > 0.95, open stars) and EROs (grey dots). 
Predicted colors for an elliptical cluster galaxy are shown 
with dashed lines, and the direction of reddening is in- 
dicated by an arrow. Evolutionary tracks are overplotted 
as in Fig. |21 to facilitate the comparison: E galaxy (red 
solid line), a low-metallicity starburst (thick solid black 
line) and a local starburst from Kinney et al. 1996 (dot- 
ted line). Error bars are typically < 0.1 mags for these 
objects. 

surements are contaminated. In the latter case, we found 
stars close to saturation in one or several filters. Thus, 
except for these understandable cases, no bright objects 
{H < 22.5) are found to fulfill our two color selection cri- 
teria. 

6.2. Catalog of optical dropouts 

Optical dropouts are defined as objects non-detected in 
all the available optical images, from the U band to the 
z band. A source is considered as non-detected in a given 
band when its magnitude corresponds to a flux below the 
Ict detection level within a 1.5" aperture. These are the 
magnitude limits reported in Tabled also used in all ta- 
bles and figures through the paper to derive limits in color. 

In the first step, we used S Extractor in "double- 
image" mode, with the H band detection images described 
above as a reference (see Sect. I5.lj) . We cross-checked the 
detection (or non-detection) of sources on their original 
images, before any geometric correction or resampling. 
Positions for all the sources detected in the H band were 
computed in the original images using the appropriate ge- 
ometrical transformations. An input catalog was created 
with these coordinates, used by SExtractor as input for 
the detection (or non-detection) of each object on the orig- 
inal image. 

The automatic procedure described above provided a 
first catalog of optical dropouts containing 122 and 38 
objects up to = 24.0, in Abell 1835 and AC114 respec- 



tively. Since we detected sources with SExtractor at very 
faint limits in flux, all the optical dropouts were carefully 
examined to reject both spurious detections in the near-IR 
bands and false non-detections in the optical bands, using 
the original images. 

Some objects in the flrst catalog were rejected because 
they were found to be located in noisy regions, close to the 
limits of the deep images or close to the haloes of bright 
galaxies. Some of them were contaminated because they 
lie too close to bright objects, in particular towards the 
cluster center. All these objects were removed from the 
dropout catalog either because they were more likely false 
detections, or because their photometry was highly con- 
taminated. A mask was created to remove the remaining 
noisy regions (bright galaxies, galaxy haloes, ...) from the 
subsequent analysis. The region masked typically corre- 
sponds to ~ 20 % of the surface in both clusters. At the 
end of the visual inspection, only ~ 20 % of the original 
sample remains in the list. 

Several dropout sources were found to be brighter in 
the Ks band image than in the H band. Their centroid de- 
termined by SExtractor on the detection H band image 
is, in principle, less accurate than in the Ks band, thus po- 
tentially leading to less optimal colors. There are 7 objects 
of this kind in Abell 1835 and 1 in AC114. We corrected 
the photometry and centroid positions for these objects by 
running SExtractor with the Ks band image as detection 
frame, keeping all other parameters unchanged. 

Since photometry was obtained in double mode, mag- 
nitudes measured by SExtractor in the infrared bands 
could be incorrect due to flux contributions at the limits 
of the aperture whereas no object is clearly seen in the 
center. For 5 objects in Abell 1835 and 2 in AC114, a 
non-detection was forced after visual inspection. 

We have considered that sources detected in at least 
two near-IR bands had more significance, since the prob- 
ability of false-positive detections in two different bands 
decreases strongly compared to our estimates only based 
on the detection band (Sect. I5.3|l . These objects consti- 
tute our "second-category" sample. Among them, we de- 
fine a "first-category" subsample including only the best- 
detected sources (having Amn < 0.4, equivalent to 2.5a 
detection within the aperture). Objects clearly detected in 
the reference H band, after visual inspection by two differ- 
ent persons, but not detected in another filter, constitute 
the "third category" sample. The remaining ones (only 
detected in H band, and dubious after visual inspection), 
are considered as a "fourth category" sample, which is 
not discussed hereafter. Tables ITT^ for Abell 1835 and l(].3l 
for AC 114 provide the coordinates and photometric prop- 
erties for all optical-dropouts from the first, second and 
third-category samples. Identification numbers increase 
with measured H band magnitudes for a given cluster. 
The number of first/second/third-category dropouts for 
Abell 1835 and AC114 is 11/7/5 and 4/4/2 respectively 
In the case of Abell 1835, the best limit for the optical 
non-detection is provided by the z band {zab ^ 27.3, cf. 
Table 0. For AC 114 the strongest non-detection criterion 
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is in R/HST with Rab > 28.-28.9, depending on the 
source location. 

Figures IC.ll and IC.2I display the thumbnail images 
of the relevant optical dropouts in Abell 1835 and AC114 
respectively. For each source, the available near-IR images 
(SZ J H Ks for AbeU 1835 and J H Ks for AC114), to- 
gether with the strongest non-detection band (z for Abell 
1835 and R/HST for AC114) are shown. 

In order to derive global properties of the different 
types of candidates, we carefully estimated our sample 
completeness and the number of false-positive detections, 
as detailed in App. ^ According to our results, false- 
positive detections are not expected up to if = 23.0, and 
they account for less than 30 % (typically 12% in A1835 
and 25-33% in AC114) for first and second-category can- 
didates brighter than H = 23.3, depending on the set of 
filters where an object is being detected (see Table IB.1|) . 

We also performed a number of additional tests on the 
reliability of optical dropouts, reported in App. [0 

6.3. Third and Fourth category candidates 

The manual classification of objects detected only in the 
reference H band into "third" and "fourth" category 
dropouts seems arbitrary at this point. We have used 
rather conservative criteria in this study to avoid the 
sample being dominated by false-positive detections, 
increasing with magnitude. However, a fraction of these 
rejected sources is actually real, although difficult to 
quantify with present data. A good example is A1835#35, 
a source kept as a secondary target for spectroscopy, for 
which we have obtained a spectroscopic confirmation 
of z = 1.68, using and [Oiii]A,A 4959,5007 lines 
detected in the J band with ISAAC (see Richard et al. 
2003 and discussion in Sect. l8.t)|) . This object should have 
been removed from our present sample: it is marginally 
detected in H and K bands, and is fainter than any other 
object reported in Tables \CT2\ and [0131 



6.4. Crude redshift estimate of high-z— candidates 
from near-IR colors 

We have produced a first list of potential "high redshift" 
candidates by applying the color-color selection criteria 
described in Sect. [3 to the catalog of optical dropouts 
(see Fig. (SJ. We select objects with a fairly red color at 
wavelengths close to the Lyman-a break/ Gunn-Peterson 
trough ("dropout"), and blue colors longward of it, indica- 
tive of a blue UV restframe spectrum. 

In practice, depending on the available near-IR pho- 
tometry, the candidates can be classified in three approxi- 
mate redshift bins between 6 and 10. From the JHK color- 
color diagram, available for both clusters, we selected a 
sample of candidates in the range 8 < z < 10. The selec- 
tion region we used is defined by : 

{H -K)< 1.0 



and [J -H)> 0.5 

and [J -H)> 1.625 {H - K) + 0.175. 

As shown in Fig. |51 the majority of optical dropouts in 
both clusters fulfill the high-z requirements. The majority 
of those in the remaining part of this diagram fulfill the 
EROs selection criterion oi R — K > 5.6, and thus they 
are possible intermediate-redshift reddened starbursts. 

To further distinguish the objects at z ^ 8 we use, 
where available, the SZJH color-color diagram to select 
candidates in the range 7 < z < 8.5 and the zSZJ dia- 
gram in the range 6 < z < 7.5. The selection in the SZJH 
diagram is defined by 

[J -H)< 2.0 
and {SZ - J) > 0.4 
and {SZ - J) > 0.8 {J - H) + 0.4 

In the zSZJ diagram it is defined by 

[SZ ~J) < 1.4 
and (z - SZ) > 0.6 
and (z - SZ) > 1.23 {SZ - J) + 0.477 

Figure El presents the color-color diagrams for the 18 
(8) first and second-category dropouts detected in Abell 
1835 (AC114). We used the location in these diagrams 
to distribute the candidates within the different redshift 
ranges. For a subsample of these optical dropouts, indi- 
vidual photometric redshifts were derived from their pho- 
tometric SEDs. The attribution of a redshift to each can- 
didate is discussed in Sect.H 

7. Results 

We shall now present the results concerning the magnifica- 
tion of the high-z candidates, their SEDs and photometric 
redshifts. Some individual objects deserve particular dis- 
cussion. We also provide some elements to understand the 
differences found between the two clusters. 

7.1. Magnification of the high-z candidates 

High-z candidates were selected based only on their pho- 
tometric properties. Their positions with respect to the 
critical lines were not considered as a selection criterion. 
However, objects located close to the high-z critical lines 
are of greater interest, because of the larger magnification. 

Figures ^1 and ^2 show the final-processed H band 
images used for the object detection with SExtractor, to- 
gether with the location of our candidates in both clusters. 
Also plotted are the critical lines at z: = 1.5 and z = 10, 
and contours of iso-magnification assuming a source red- 
shift of z = 9. computed from the lensing models for Abell 
1835 (similar to Smith et al. 2005) and AC114 (Natarajan 
et al. 1998, Campusano et al. 2001). The position of these 
lines is weakly sensitive to source redshift within the rele- 
vant range z ~ 6 to 10. For a given redshift estimate, the 
location of the high-z candidates on the field allows us 
to derive their magnification factors (see Tables and 
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Fig. 9. Color-color diagrams showing the location of all optical-dropouts detected in Abell 1835 (a-b-c), AC114 (d), 
and the delimitation of the selection region used in the different redshift domains: JHKs (z ~ 8 — 10), SZJH 
{z ^ 7 ~ 8.5), and zSZJ (z ~ 6 — 7.5) (see text). The identification numbers are indicated according to Tables Rr2l 
and l(].3l For a given diagram, circles and squares correspond to objects detected in three and two filters respectively; 
upper limits are displayed by an arrow. Optical dropouts fulfilling the EROs definition are shown in grey. Dropouts 
are all detected in the H band, and non-detected in the z band. Sources are presented in diagrams (a) and (d) if they 
are detected either in J or K, and in diagrams (b) and (c) if detected in J or SZ. 



IC.3II . /ig and /iio give the magnification factors assuming 
a source redshift of 6 and 10 respectively. Where appli- 
cable, the adopted magnification // is computed assuming 
the "adopted" redshift z described in Sect. 17.21 



tiple image system. Because of the relative insensitivity 
to source redshift and position on the image plane, a re- 
fined version of lensing models will not change our present 
results and conclusions. 



The uncertainty in the magnification factor associ- 
ated with the uncertainty on the redshift value is usually 
smaller than 10%, except for a few objects exhibiting the 
largest magnification factors (/i > 10), i.e. located within 
a few arcsecs of the critical lines. Also the magnification 
factor at a given position on the image plane varies slowly 
with redshift for sources located more than ^ 10" away 
from the critical lines. Since lensing models are mainly 
based on the identification of multiple images with secure 
spectroscopic / photometric redshifts, the uncertainty in 
the derived magnification factor is usually smaller close to 
these regions. However, it could be larger in the case of 
Abell 1835, because the model is based only on one mul- 



For objects located close to the critical lines at high 
redshift, we used the same lensing models to look for pos- 
sible multiple images which could affect our analysis (num- 
ber counts, etc.) or allow us to better constrain the posi- 
tion of the critical lines at z > 6. With the present data, 
we could not find any pair of objects that would be a fair 
multiple image candidate, although this possibility cannot 
be ruled out because of completeness considerations (see 
Sect. ESI . 
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Fig. 10. Left : H band image of the lensing cluster Abell 1835 showing the location of the critical lines at z = 1.5 
(thin solid curve) and z — 10 (thick solid curve). All candidates are shown with crosses; identification numbers are the 
same as in Table lU^ First/Second-category dropouts are circled and fourth-category objects are not labeled. Right: 
location of the same objects relative to the magnification across the field. Contours are overplotted for magnification 
values of 1, 2 and 3 magnitudes, computed assuming sources at z = 9, although the position of these lines is weakly 
sensitive to source redshift within the relevant range z 6 to 10. White lines delimit the footprint of the i?702 WFPC 
image, covering ^ 46% of the whole ISAAC field of view. 




Fig. 11. Same caption as Fig. ^|for the cluster AC114. Identification numbers are the same as in Table 1^751 White 
lines delimit the footprint of the i?702 WFPC image, covering ^ 85% of the whole ISAAC field of view. 



7.2. SED properties and photometric redshift estimates 

The position of optical dropouts on the different color- 
color diagrams provides a first estimate of their photo- 
metric redshifts, and an objective criterion allowing us to 
classify them into different redshift intervals. This crite- 
rion can be refined for about 30% of our candidates pre- 
sented below, for which the S/N is sufhcient to derive indi- 
vidual probability distributions in redshift. These sources 
are among the brightest candidates in the H band, or 
have been detected in other filters with S/N~ 4 — 10. 



Of particular interest are the optical dropouts which can 
be unambiguously excluded from the z > 6 sample using 
photometric redshift considerations. 

Probability distributions and photometric redshifts 
were derived for our candidates from broad-band photom- 
etry over a wide wavelength interval, using an adapted ver- 
sion of the public photometric redshift software Hyperz 
(Bolzonella et al. 2000). Best-fit redshifts and redshift 
probability distributions between z = and 12 were 
computed through a standard SED fitting procedure. We 
used a variety of template models: starbursts (Kinney et 
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al. 1996; SBS0335-052, Izotov 2001), evolutionary mod- 
els from the GISSEL library (Bruzual & Chariot 1993), 
empirical templates from Coleman et al. (1980), and the- 
oretical templates for zero metallicity (PopIII) and low 
metallicity starbursts (Schaerer 2003). Intrinsic redden- 
ing was considered as a free parameter ranging between 
Ay = and 3 magnitudes, according to the Calzetti et 
al. (2000) extinction law. The Lyman forest blanketing is 
included following the prescription of Madau (1995). The 
non-detection in the optical bands was used as a constraint 
when computing photometric redshifts. 

For further discussion on the candidates, we shall 
adopt a redshift value z combining all the constraints ob- 
tained from the color-color diagrams and the photometric 
redshift determinations. Redshift values z and redshift in- 
tervals adopted for each optical dropout are reported in 
Tables IC.2l and IC.3I Depending on its photometric red- 
shift probability distribution P{z), we attribute a redshift 
quality to each candidate (0z in Tables IC.2I and IC.3|) as 
follows: 

(A) Objects displaying a unique solution in their probabil- 
ity distribution in redshift (P(z)), irrespective of the 
redshift, with a good SED fitting for this unique solu- 
tion (i.e. absolute probabihty higher than 90% in most 
cases) , with no other secondary solutions with reduced 

better than x^(best fit) -t- 1. The best fit redshift in 
this case corresponds to the 5, and zl,z2 given in the 
tables correspond to Icr redshift intervals (68% confi- 
dence intervals). 

(B) Objects showing a degenerate solution between z > 
6 — 12 and a lower redshift solution, in general within 
z ^ 1.5-2.5. In this case, the two solutions are equally 
significant, with a good SED fitting in both cases (re- 
duced ^ !)• Wc adopt in this case the higher red- 
shift solution in a similar way as for (A). 

(C) Objects for which no reliable individual photomet- 
ric redshift solution could be obtained, either because 
their P{z) is basically "flat" , without a significant solu- 
tion, or because multiple degenerate solutions exist, all 
of them providing a poor fit in terms of absolute prob- 
ability. In this case, the redshift limits [zl, z2] and z 
in the Tables correspond to the mean redshift and the 
redshift intervals defined by the color-color diagram 
selection: 

[6.0-7.5] and 6.75 for zSZJ 
[7.0-8.5] and 7.75 for SZJH 
[8.0-10.0] and 9.0 for JHK 

For objects displaying a peak at high-z in their probability 
distributions P{z) (quality types A and B), the best-fit z 
is always found as expected within the redshift interval 
defined by the color-color diagrams. SZ band photometry 
is not available for AC114, and thus we only considered the 
[6.0-8.0] and [8.0-10.0] redshift intervals for this cluster. 



Figures Kill and KL 21 also display the SEDs and best- 
fit models for optical dropouts, only for first and well- 
detected second-category sources. Several of these objects 
seem too bright to be high-z sources. We discuss below 
these and other peculiar objects individually. A more de- 
tailed description of individual objects will be presented 
elsewhere. 

7.3. Contamination by mid-z interlopers 

Because of the large photometric errors of the selected op- 
tical dropouts in the near-infrared bands (up to ^ 0.6 mag- 
nitudes in some cases) our high redshift sample is suscep- 
tible to contamination by low redshift interlopers falling 
in our color-color selection regions. We derive a basic idea 
of this contamination factor based on the spectroscopic 
catalog from the HDFS (Noll et al. 2004) for which the 
near-infrared photometry was obtained with the same in- 
strument and filters (Labbe et al. 2003). 

Photometric errors were introduced in the HDFS cata- 
log following a Gaussian distribution of fixed a in each J, 
H and Kg band. Different values of a were used to mimic 
the typical S/N in the sample of optical dropouts, and 
different redshift intervals were considered for the HDFS 
sample. The worst contamination level by mid-z interlop- 
ers on the [J—H) vs {H~K) color-color selection diagram 
ranges between 20 and 25%, depending on the redshift 
intervals, when applying a 0.6 mag. photometric scatter 
to all filters, i.e. error bars sensibly larger than those re- 
ported in Tables IC.2I and IC.3I None of the contaminant 
sources would have been selected as an optical dropout. 
Therefore, this source of contamination should represent 
a second order correction for our sample. 

7.4. Individual objects in Abell 1835 

Optical dropouts in this cluster are distributed as fol- 
lows (cf. Table [U^l Figure EH). We have detected 7 
objects satisfying the ERO criterion. Among them, two 
(A1835#2 and A1835#17) are unambiguously identified 
as low-z sources, two display a non-standard behavior in 
then SED (A1835#10 and A1835#11), and one is clearly 
variable in (at least) the SZ band (A1835#4). In general, 
optical dropouts displaying anomalous SEDs as compared 
to young starbursts are either variable sources (A1835#4) 
or sources whose nature could not be determined with the 
present data (A1835#8, A1835#10, A1835#ll). These 
objects, marked as "Ex." in Tables Rr2l and IC.3I are not 
considered as true high-z candidates and are therefore ex- 
cluded from the high-z sample. Similar arguments can be 
used to remove the other two "bright" EROs (A1835#1, 
A1835#3) from the sample. 

A general comment concerning the "bright" EROs is 
needed here. As mentioned above, only 2 of these in A1835 
and 1 in AC 114 could be identified unambiguously as 
intermediate-z galaxies using photometric redshift con- 
siderations. All the others are difficult to reconcile with 
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normal "mid-z" galaxies because of their relatively blue 
continuum in the near-IR JHK domain, yielding solu- 
tions at z ~ 6-8 which are equally likely or even better 
(see Fig. IC.ll and Table lH^ . However, they are extremely 
luminous if at such high-z, typically Mb ~ —25.0 to —26.5 
for the brightest ones. Although "bright" high-z sources 
of this kind may exist, up to 50L*, as recently proposed 
by Mobasher et al. (2005) for their massive post-starburst 
galaxy at z~6.5, we exclude these sources from the high-z 
sample, and consider the low-z solution as more plausible 
at this stage. Spectroscopic information is needed to fur- 
ther characterise these faint lensed EROs, and determine 
their nature by a measurement of their redshifts. 
Some sources in A1835 deserve specific comments. 

— A1835#2 : This source corresponds to J5, the 
near-IR counterpart of the SCUBA-selected galaxy 
SMMJ14009-I-0252 (Ivison et al. 2000, Small et al. 
2002). The near-IR photometry reported here is com- 
patible with the recent results by Frayer et al. (2004). 
It satisfies the ERO criterion, and it is likely a low-z 
source as suggested by our SED fitting result: z=1.34 
(with la ranging between z=1.18 and 1.64). This solu- 
tion is in good agreement with the redshift constraints 
O-"^ ^ -2 ^ 1-8 derived from the radio-submm spectral 
index af^4 , but inconsistent with the range 3 < z < 5 
suggested by the submm colors (Ivison et al. 2000). 

— A1835#4 : As mentioned in Appendix O this source 
displays a difference of 0.4 mags in SZ (more than 3 
a) between the two series of images obtained on the 
19 April and the 15 May, which seems to indicate an 
intrinsically variable source. 

— A1835#8 : This object (previously named 
A1835#1916, also known as the z ~ 10 candidate) 
was studied in detail in Pello et al. (2004ab). The 
photometry presented here is an improved version; 
the present and earlier magnitudes in the common 
filters are compatible within la error-bars. The field 
around A1835#1916 has been reobserved between 
30 May and 6 June 2004 by Bremer et al. (2004) 
with NIRI/GEMINI in the H band. Surprisingly, the 
object is not redetected in these images, which are at 
least ~ 0.5 mag deeper than the ISAAC images taken 
approx. 15 months earlier. The reality of our initial 
photometric detections is not questioned by Bremer 
et al. who reconfirm it using our data, although the 
photometric properties of this source are still a matter 
of debate (Lehnert et al. 2005, Smith et al. 2006). The 
detection in 3 bands where the object is re-detected, 
including our new SZ images (see Table ICT^ . makes a 
spurious event highly unlikely (about 12% probability, 
from our estimates given in Table IB.ip . When we 
consider the present results together with our previous 
findings (the source was virtually non-detected in 
our J images), and the non-detection by Bremer et 
al. in the H band with new independent data, this 
source could be intrinsically variable. Its nature (and 
hence also its redshift) presents a puzzle, and will be 



discussed elsewhere. Hereafter we do not consider this 
source within the high-z sample. 

— A1835#35 : This is a third-category candidate (previ- 
ously named A1835#2582) and a rather unusual emis- 
sion line galaxy already studied in detail by Richard 
et al. (2003). As for #8, the present and earlier mag- 
nitudes in the common filters are compatible within 
la error-bars. The marginal detection of this object in 
the / and J bands reported by Richard et al. (2003) is 
found to be non-significant with our new and more con- 
servative error estimates. Thanks to [Oiii] AA4959,5007 
and H/3 detected in the J band, #35 has been iden- 
tified as an extremely faint galaxy at z = 1.68, with 
Mb ^ —16.4 and a gravitational magnification of ~ 2 
magnitudes. This object has been removed from the 
photometric sample of high-z candidates. 

7.5. Individual objects in AC114 

From the 10 optical dropouts selected in AC114 (cf. Table 
IC.3I Figure [(T2|l . 8 are in the first and second-category 
samples. Among those, only one object satisfies the ERO 
criterion (AC114#1, cf. below). The lack of z and SZ 
photometric data for this cluster precludes a further clas- 
sification into redshift bins between z ~ 6 — 8. 

— AC114#1 : This source satisfies the ERO criterion, and 
it is likely a low-z galaxy according to our SED fitting 
result: z = 1.62 (with la ranging between z — 1.58 
and 1.89). 

7.6. Differences between the two lensing fields 

The total number of remaining first and second-category 
high-z candidates in AC114 as compared to Abell 1835 is 
found to be in a ratio of 7/10 (9/15 whith third-category 
candidates), after excluding EROs in both two clusters. 
Several reasons could explain this difference, in addition 
to field-to-field variance: 

— The depth of the near-IR images is less in AC 114 as 
compared to Abell 1835, by 0.1 and 0.4 magnitudes 
in J and Ks respectively. Also the overall detection 
image in H is noisier for AC 114, as shown in Table 
IB. II rendering the identification and visual inspection 
of near-IR detections more difficult. We expect a larger 
fraction of (blind) false positive detections in AC114 
than in Abell 1835 f Table IrT|| . Thus, if the sample 
was dominated by such detections, the number of can- 
didates should be smaller in Abell 1835 than in AC114. 
However, the opposite trend is observed, which means 
an efficient (manual) control of the sample. 

— The optical images used to identify dropouts are not 
identical in the two cases. In Abell 1835, the main 
constraint comes from a deep z band image (limiting 
magnitude zab ^ 27.3, whereas it is Tab ^ 27.2 in 
AC114). In AC114, the strongest non-detection crite- 
rion is set by the R/HST image (limiting magnitude 
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Rab ^ 28.7 — 28.9 in the relevant region of the field, 
whereas it is Rab ^ 28.0 in Abell 1835). Star-forming 
and highly-reddened intermediate-z sources such as 
A1835#35 could survive more easily in the Abell 1835 
sample than in AC 114 because of the difference in the 
R band. 

— Geometrical considerations coming from lensing are 
also to be taken into account. Abell 1835 is not cen- 
tered on the CD galaxy, and the two clusters have dif- 
ferent redshifts (z = 0.252 and 0.312). As shown in 
Fig. [TUIand [TTl the magnification factors across the 
ISAAC field of view are different thus leading to differ- 
ent effective surveys, as discussed in detail in Sect.|Hl In 
summary, within the redshift domain considered here 
(6 ^ ^ ^ 10), the two clusters cover about the same 
area on the sky close to the critical lines for the largest 
magnifications factors beyond 5 (24 and 21% of the 
total surface respectively), and thus we expect (and 
observe) the same number of such sources in the two 
fields. On the contrary, for magnification factors below 
5, the two clusters exhibit a different behavior: 31 and 
47% of the total surface respectively for Abell 1835 
and AC 114, for magnifications factors between 2 and 
5, and 45 and 32% of the total surface respectively for 
magnifications factors below 2 (0.75 magnitudes). In 
other words, the sample observed in Abell 1835 is ex- 
pected to be dominated by sources with magnification 
factors below 2 (up to ^50-60%, assuming sources with 
a uniform distribution in z and same luminosities, af- 
ter completeness correction), whereas the majority of 
candidates detected in AC114 are expected to exhibit 
magnification factors above 2 under the same condi- 
tions. This is the qualitative behavior of candidates 
actually observed. For magnifications factors above 2, 
and taking into account that we are dealing with small 
number statistics, the relative number of candidates 
detected in both clusters is in good agreement with ex- 
pectations. On the contrary, the simple considerations 
given here cannot explain the excess in the number 
of candidates detected in Abell 1835 with respect to 
AC114, for magnification factors below 2. 

8. Discussion 

In this section we discuss the implications of the present 
results in terms of abundance of star-forming galaxies at 
6^-2^ 10. We present and discuss the intrinsic properties 
of candidates actually detected, after correction for lens- 
ing magnification, and we compare the observed counts 
with order-of-magnitude expectations obtained from sim- 
ple modeling. The observed number densities of candi- 
dates, as seen through gravitational lenses, are translated 
into effective number densities through a careful modeling 
of lensing effects, easily comparable with blank field stud- 
ies. The luminosity function and the cosmic star formation 
rate derived from our sample of 6 < z < 10 candidates is 
also presented and discussed. The large correction factors 
applied to this sample make the determination of inte- 



grated quantities, such as luminosity functions and SFR 
densities extremely difficult. In the final section, we briefly 
describe the preliminary results obtained on the spectro- 
scopic follow-up of the photometric candidates. 

8.1. Intrinsic properties 

The typical magnification values of our candidates range 
between 1.5 (~ 0.44 mags) and 10 (2.5 mags). For some 
objects very close to the critical lines, we found magnifi- 
cation values n > 25. However, because of the underlying 
error in the precise location of the critical lines from the 
models, we prefer to adopt a more conservative lower limit 
of /z = 25 for these objects. The average(median) magni- 
fication values among the first-priority high-z candidates 
are 6.5 (2.3) in AbeU 1835 and 7.9 (3.5) in AC114. 

We derived the unlensed iisoo luminosity, at 1500 
A in the restframe, for all high-z candidates, using the 
adopted value z for the redshift and a fiux estimate from 
the photometry in the band closest to this restframe wave- 
length {SZ, J or H). Deriving the intrinsic properties of 
6 < z < 10 galaxies can be challenging due to the signifi- 
cant level of contamination. 

After correction by the lensing magnification affecting 
each object (see Sect. 17.111 the unlensed iisoo luminosities 
were converted into Star Formation Rate (SFR) through 
the usual calibration from Kcnnicutt (1998): 

SFR (Mq yr-i) = 1.05 lO'^" L1500 (ergs s'^ A"') 

These physical properties are summarized for each can- 
didate in TableESlandinil The typical SFR obtained for 
objects included in the final sample (excluding EROs and 
anomalous SEDs) is about ~ 10 M© yr"^, with extreme 
values ranging between a few units and ^ 20 M© yr~^. 
The conversion of iisoo into SFR assumes a constant star 
formation at equilibrium, and such conditions arc not nec- 
essarily reached in these objects. 

Interestingly, although the selection criteria are only 
based on near-IR colors irrespective of magnitude, almost 
all the photometric candidates fulfilling our selection crite- 
ria are fainter than H — 23.0 (AB ^ 24.5). Only three ex- 
ceptions are found in Abell 1835 among the possible low-z 
EROs, as described above. After correction for magnifica- 
tion accross these fields, the lack of "bright" sources means 
that we have not detected young starbursts at z ^ 6 — 10 
more luminous than typically L1500 = 3 x 10"'*^ ergs s~^, 
i.e. more massive than typically a few 10* M© (for star- 
bursts younger than 10^ yr, under standard assumptions 
for the IMF). 

Also, a direct comparison between low-z and high-z 
SEDs shows that our high-z candidates tend to be very 
blue in H — K. For z > 6 candidates, H — K colors pro- 
vide an estimate of the restframe UV slope /?, assuming 
3ct detection levels in K when the source is not detected 
in this band, or when the S/N in this band is lower than 
3(7, with large uncertainties due to photometric error bars. 
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The UV slope f3 usually ranges from —2.5 to —3.5 for the 
first category dropouts, with two sources (A1835#7 and 
AC114#2) reaching ~ —3.9 at la, and between -0.7 and 
—3.0 for the second category dropouts. This systematic 
trend towards extremely blue colors was also reported by 
Bouwens et al. (2004b) for their sample of z ~ 7 — 8 can- 
didates. 

Although the optical dropouts we found are stretched 
by the magnification factor fi, they appear as point-like 
sources in our ground-based images. If we assume a mini- 
mum magnification of 0.4 mags for all the field, the phys- 
ical size of these objects at z > 7 is likely to be smaller 
than 1.7 kpc, in good agreement with size calculations of 
Kneib et al. (2004) and Bouwens et al. (2004b). 

8.2. Observed number counts vs model expectations 

The efficiency of using lensing clusters as gravitational 
telescopes to find high-z galaxies can be evaluated with 
model expectations and simple assumptions. 

We first use a simple model to estimate the expected 
number counts, both in blank fields and lensing clusters, 
taking into account our photometric selection criteria. We 
use semi-analytical models for dark-matter halo formation 
(the Press-Schechter formalism. Press & Schechter 1974). 
Starburst models presented in Sect. [S] were used to scale 
the SED, assuming that the fraction of the baryonic mass 
converted into stars is fixed to /, =0.1 within the redshift 
interval 6 < z < 11. We focus on two different "extreme" 
assumptions for the IMF: a standard Salpeter IMF from 1 
to 100 M0 and a top-heavy IMF (with stars ranging from 
50 to 500 Mq). 

To correct these estimations for the relative visibility 
time of the bursts, which are typically of t* = 10^ yr 
timescales from our simulations, we multiply the number 
density of objects expected at a given redshift z by the 
duty-cycle factor : t^ (H-z)/(t// (z) — t^f (17)), where tniz) 
is the age of the Universe at redshift z. This corresponds 
to the probability for a burst to be visible at redshift z, as- 
suming that all haloes convert a constant fraction of their 
baryonic mass into stars, at some time between redshifts 
17 and z. 

Lensing introduces two opposite trends in the observed 
sample compared to blank fields: 1) gravitational magnifi- 
cation by a factor /i, increasing the number of faint sources 
and thus the total number of sources, and 2) reduction 
of the effective surface by the same factor thus leading 
to a dilution in observed counts. We explicitly compute 
the expected number counts with these models in clus- 
ters by a pixel-to-pixel integration of the magnification 
maps as a function of redshift, using the lensing models, 
after masking all the pixels lying in the mask of bright ob- 
jects described in Sect. 16.21 The expected number counts 
up to if < 24 are very similar for the two clusters (less 
than 10% difference over the redshift interval), and thus 
we averaged together both results into a unique "strong 
lensing cluster" prediction. These results hold for interme- 



diate redshift clusters (z ~ 0.2 — 0.3) and should be revised 
for lensing clusters at higher redshifts and/or a different 
field of view, although the trends remain the same. 

The comparison between expected and observed num- 
ber counts of galaxies in the field of ISAAC, up to < 24, 
per redshift bin Az — 1, in a blank field and in the 
field of a strong lensing cluster are presented in Fig. 1121 
Blank field number counts are not corrected for bright- 
object masking, and thus they correspond to an upper 
limit. Observed number counts in the two lensing clus- 
ters have been corrected for photometric incompleteness 
(typically a factor of 0.05 — 0.1) using the simulations re- 
ported in Sect. 15.31 sample incompleteness (by a constant 
factor Csampie ~ 0.42 as detailed in Appendix 0) and 
for the expected fraction of false-positive detections, from 
our estimates given in Table IBJI Results are presented in 

Fig. [121 

For each redshift bin, we corrected the number counts 
for both first and second-priority candidates (using z as 
their redshift estimate), and the results are directly plot- 
ted on the model expectations, without any normalization. 
Observed number counts are upper limits, as our sample 
is likely contaminated by low-z interlopers. 

As shown in Fig. ^1 strong lensing fields are a factor 
of 5 — 10 more efficient than blank fields of the same size 
in the z ~ 7 — 11 domain, all the other conditions being 
the same. Observed number counts of candidates at z '--^ 
6 — 8 are in good agreement with these order of magnitude 
estimates, in particular with ^10% of the baryonic mass 
converted into stars at z > 6. At z > 8, the observed 
number counts are more consistent with a top-heavy IMF, 
with a large cluster-to-cluster variance. 

8.3. Lens-corrected number densities of high-z sources 

We have used the lensing models to derive the effective 
areas and corresponding volumes surveyed in the different 
source planes. The aim is to translate the observed num- 
ber densities of candidates into effective number densities 
easily comparable with blank field studies. We also correct 
our observed sample of candidates for incompleteness and 
false-postive detections. 

Magnification and dilution effects by the lensing field 
are carefully taken into account to compute number den- 
sities and derived quantities. The average magnification 
value over a whole ISAAC field is about 2, thus leading to 
a dilution close to 50% over the whole field. However, a 
careful modeling is needed to properly take into account 
the intrinsic incompleteness of the sample as a function of 
redshift and position on the field. For each candidate in the 
field, with observed magnitude Ho, we compute the mag- 
nification factor A/(f2, z) as a function of its position and 
redshift z, as well as the lens-corrected magnitude (re- 
ferred hereafter as effective magnitude), using the lensing 
models presented in Sect. 17.11 The effective completeness 
T] {He,z) gives the ratio between the observed number 
counts in the lensing field, No{He,z), and the equivalent 
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Fig. 12. Comparison between the expected number 
counts of galaxies in the field of ISAAC, up to H <24, 
per redshift bin Az = 1, in a blank field and in the field 
of a strong lensing cluster (see text for details) . Expected 
counts are obtained with the simple model discussed in 
Sect. 18.21 for two extreme IMF: a standard Salpeter IMF 
(lower curves) and a top-heavy IMF (upper curves). The 
differences between a blank-field (dashed line) and lensing 
fields (solid lines) are more pronounced at higher redshift. 
Observed counts are displayed for the two lensing clusters, 
corrected for incompleteness effects. Solid lines display the 
results for first-category sources only, whereas dotted lines 
correspond to first and second-category candidates. 

value measured on a blank field of the same observed sur- 
face, N{He,z), complete up to the magnitude He'. 
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This quantity can be written as a function of the solid 
angle surveyed on the sky, A 17, 
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where C{Ho) stands for the photometric incomplete- 
ness correction for an observed H band magnitude Ho 
(plotted in Fig. ^ and Csampie — 0.42 is the addi- 
tional incompleteness factor for the sample, as given in 
Appendix IbI 

In practice, we computed the values of 77 as a func- 
tion of _ffe and redshift using the magnification maps over 
the field, after masking all the pixels lying in the mask 
of bright objects described in Sect. 16.21 Figure IT^ shows, 
for each cluster, the location of the candidates in the 
{z,He) plane. Overplotted are different models of high- 
z starbursts, and completeness levels ry. Excluding the 3 
brightest EROs in Abell 1835 (which are possibly low-z 



sources), all our candidates range from He— 24.0 to 27.2 
(AB ~ 25.5 to 28.7). From these diagrams, it appears that 
our sample of candidates is sensitive to stellar mass scales 
in the range 10^ — 10® M©, and that our typical correction 
levels range from 1 to 15% (including both lensing dilu- 
tion and photometric incompleteness). Candidates with 
the smallest t] factors in these diagrams have the largest 
weights in the number densities and derived quantities. 

8.4. Luminosity Function 

A fair estimate of the luminosity function (LF) at 1500 A 
can be derived from the intrinsic luminosities of our can- 
didates presented in Sect. 18.11 using the same approach 
as in the previous section to compute number densities 
corrected for incompleteness and spurious detections. As 
discussed in the previous sections, the correction factors 
applied to this sample are relatively large, thus leading to 
large error bars in the LF determination. We discuss here 
two redshift intervals for which we have enough sources 
for this exercice: 6 < 2; < 10 and 8 < 2; < 10. In practice, 
we derive the number density of objects in the co- volume 
surveyed, with Alogio{L) = 1, after correction for the in- 
dividual 77 factors. Error bars are estimated as a combina- 
tion of two independent sources of noise: the la confidence 
levels for a Poisson distribution and the uncertainty in lu- 
minosity introduced by the rj factor (typically a factor of 
3). The latter is a combination of the photometric error 
bars, and uncertainties in the incompleteness corrections 
and lensing modeling. When no object was detected in a 
luminosity interval, we corrected the Poisson-noise upper 
limit by the typical effective completeness rj for this lumi- 
nosity. We also corrected our data points for the fraction 
of false-positive detections expected from Table IB. II 

The combined -L1500 Luminosity Functions for both 
clusters, with the corresponding error-bars, are given in 
Fig. Only first and second-priority candidates have 
been considered, but the difference obtained when using 
the full sample is within the la error bars. 

The observed LFs have been fitted by the STY method 
(Sandage, Tammann & Yahil 1979), a maximum likelihood 
fit of data points to the analytical Schechter function <f){L) 
(Schechter et al. 1976): 



</.(L)di = 0^500 — 



exp 



-- 4- 

L* \L* 



assuming that this function provides a good representa- 
tion of the data. Due to a lack of information towards the 
faint end, a strong degeneracy is expected between L* and 
a, which we do not discuss here. To avoid this problem, 
we assumed a fixed value of a = 1.6, corresponding to 
the Steidel et al. (1999) determination for Lyman Break 
Galaxies (LBGs) at z ~ 4, and leave the other parame- 
ters free. The typical values found for L* are 10^^'^ ergs 
s^^ s^^ A and these results are not affected by the 
way we binned the data points. The STY fits to the data 
presented in Fig. 
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Fig. 13. Effective magnitudes H^, corrected for lensing, as a function of z for the two lensing clusters. Values of 
effective completeness ri{H(,, z) with respect to a blank field (for number counts complete up to He) are overplotted 
as dotted lines for the 0.5 to 40% levels (see text for details). First and second-category (red) and third-category 
(green) candidates are positioned at their adopted redshift z, with errors bars in z corresponding to [zl — z2] range 
(see text for details). Redshift intervals considered for third-category candidates are given at the bottom of the figure. 
Photometric errors in H,, coincide with photometric errors observed in H. Solid lines display the predicted magnitudes 
versus redshift for 4 different starbursts models used in Sect. El with a stellar masses scaling to 10^ M0 (thick line) 
and 10^ M0 (thin line). From top to bottom, these models correspond to single bursts with Salpeter IMF and stellar 
masses ranging from (1) 1 and 100 Mq, and (2) 50 and 500 M0. For comparison, the dashed line corresponds to a 
constant star-forming model with age 10* yr. 



Without the blind correction for false-positive detec- 
tions, the data points increase by ~0.4-0.6 dex, which is 
a relatively small effect compared to the error bars. 

For comparison, we overplot in Fig.^lthe LF fit found 
by Steidel et al. for LBGs at z 3, after correction for 
differences in the respective cosmological parameters. It 
is shown as a thick dashed line in Fig. without any 
renormalization to fit the data points. This LF for LBGs 
at z ~ 3 seems to be slightly higher (by ~ 0.5 dex) than 
the LF for our candidates, but compatible within the la 
error bars. Our results are also fully compatible with the 
LF derived by Bouwens et al. (2005) (presented as a dotted 
line in Fig. 114(1 for their sample of z 6 candidates in the 
UDF, UDF-Ps and GOODS fields in the low-luminosity 
regime, i.e. for L1500 < 0.3L*^3, but we do not see the 
turnover observed by these authors towards the bright end 
relative to the z ^ 3 LF. 

We also display in Fig. ^1 the LF derived from the 
simple models presented in Sect. 18.21 in the redshift inter- 
vals 6 < z < 10 (red lines) and 8 < z < 10 (black lines), 
for the standard Salpeter IMF (thin dashed lines) and the 



top-heavy IMF (thin solid lines). A better overall fit to 
data points is obtained with the top-heavy IMF. 

Considering only first-priority or first-|-second priority 
candidates does not change the results substantially. Also, 
including or removing the brightest EROs in Abell 1835 
does not change these conclusions. 

8.5. Cosmic Star Formation Rate at z > 6 

We use the individual properties computed in Sect. IS.ll to 
derive the Cosmic SFR value for each redshift bin consid- 
ered in our previous analysis. The large correction factors 
applied to this sample, which dominate the error bars on 
LF measurements, make the determination of SFR densi- 
ties challenging. We use different approaches here to de- 
rive an estimate for this important quantity. On the other 
hand, since these objects are only photometric candidates, 
the obtained values are to be considered as upper limits. 
However, all values derived here neglect possible extinc- 
tion corrections. 

The first estimate is obtained in a very simple way. 
We compute the total SFR in a redshift bin by summing 
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Fig. 14. Combined I/1500 LFs for the two fields, for two 
redshift intervals: 6 < z < 10 (red) and 8 < z < 10 
(black). Data points are corrected for spurious sources, 
and error bars combine Poisson noise statistics and uncer- 
tainty in the effective completeness. Note the large correc- 
tion factors applied to this sample, which translate into 
large error bars. The STY fits to the LF data are rep- 
resented by thick solid lines. For comparison, the LF fit 
by Steidel et al. for LBGs at z ~ 3 is also overplotted 
(thick dashed line), as well as the z ~ 6 fit from Bouwens 
et al. (2005) (thick dotted line), without any additional 
renormalization to fit the data. Also the ^1590 LF derived 
from the simple models presented in Sect. is shown, 
in the redshift intervals 6 < z < 10 (red thin lines) and 
8 ^ 2 < 10 (black thin lines), for the classical salpeter 
IMF (dashed lines) and the top-heavy IMF (solid lines). 
The cosmic SFR value has been derived by integrating 
this LF down to 0.3L*^3 (shown as a vertical grey line). 
See text for more details. 

all individual contributions to the SFR within this bin, 
after correcting each object by its 77 value and its expected 
probability of being a false-positive detection (Table IB. 1(1 . 
We divide this result by the total covolume surveyed in 
this redshift bin accross the ISAAC field area, assuming 
a blank field (because magnification/dilution effects are 
already included in the rj factor). This corresponds to 
2.9 10"^ Mpc^ for z G [6-7.5], 2.6 10* Mpc^ for z G [7-8.5], 
3.8 10'' Mpc3 for z e [6-8], 3.1 10"* Mpc^ for z G [8-10], 
for the size of an ISAAC field. The resulting SFR densities 

obtained are relatively high, even when the sample is 
restricted to first and second-category candidates: p+ = 
3.31 10-3 Mq yr-i Mpc^^ for z e [6-7.5], 2.46 10"^ M© 
yr-i Mpc-3 foj.^ ^ [7-8.5], and 1.20 lO^^ for z e [8-10]. 

Another estimate of the cosmic SFR, allowing us to 
compare results to previous findings, can be obtained by 
integrating the Luminosity Function fit found in Sect. 18.41 
down to 0.3 i*=3, the same limit used by Bouwens et al. 
(2004b). In this case, we use the same redshift bins defined 
in Sect. 18. 4I and we obtain a somewhat lower value of the 
SFR for z in [6-10] if we consider the first and second 
category candidates. Considering only the first category 



candidates, these values are lower by a factor ~ 3, as can 
be seen from Table summarising our different cosmic 
SFR density estimates. 

The final results are shown in Fig. 1151 applying the cor- 
rection for false-positive detections, for comparison with 
the cosmic SFR obtained in other surveys carried out on 
blank fields, without applying any extinction correction. 
Our error bars were computed using Poisson noise statis- 
tics in the number of objects within each redshift bin. The 
numerical values of the cosmic SFR are also summarized 
in Table El 

When considering only the first-category candidates 
(i.e with a la error AniH < 0.4), our results in the 
z ~ 6 — 10 domain are compatible with previous find- 
ings. However, our estimate of the comoving SFR den- 
sity at z 8 — 10 seems to be larger than all values 
derived at z ~ 4 — 6, although compatible within the 
error-bars. Taken at face value, our findings seem to be in 
good agreement with some theoretical cosmic SFR den- 
sity models previously published; e.g. with the model of 
Barkana & Loeb (2001, their Fig. 29) for a reionization 
redshift between 6 and 8, recent hydrodynamical models 
of Nagamine et al. (2005), and with the self-consistent 
reionization models of Choudhury & Ferrara (2005). 

However, compared to recent studies in the Hubble 
UDF at similar redshifts (Bouwens et al. 2004b, 2005a) our 
SFR density, or upper limits thereof, is larger by roughly 
1 dex. This difference is related to the bright end of the 
LF, i.e. £1500 ^ 0.3L*^3. In all cases the sources are pho- 
tometric candidates therefore providing upper limits to 
the actual UV fiux densities. The effective fields surveyed 
are small in all cases, thus leading to a strong field-to- 
field variation of ~ 20 — 30% in the number of sources. 
Cluster-to-cluster variations already discussed in Sect. 17.61 
are clearly seen in our sample, although lensing and pho- 
tometric considerations could account for most of them. 
Recent spectroscopic results obtained by Le Fevre et al. 
(2005) on the VVDS, for an I-band fiux-limited sample 
of galaxies up to z 5, indicate that the Universe con- 
tains more star-forming galaxies towards the bright end of 
the LF than previously reported using color-color selection 
techniques, suggesting an active star formation activity in 
the redshift domain covered by the present survey. 

There are several ways to reconcile our SFR measure- 
ments with Bouwens et al.'s values, in addition to field- 
to-field variance. Residual contamination by false sources 
combined with lensing effects are able to affect our results 
in different ways. 

On the one hand, residual spurious sources constitute 
a potential source of contamination if the adopted correc- 
tions are underestimated. The extremely blue H-K col- 
ors obtained for the stacked images argue for a signifi- 
cant residual contamination for the faint third category 
dropouts, as discussed in detail in Appendix ICl For this 
reason, third category dropouts were not used to derive 
LFs and SFR densities. The blind corrections applied to 
the first and second-category samples to obtain the above 
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results are in good agreement with the contamination lev- 
els derived from the stacked images and the UV slope. 

On the other hand, observed number counts could be 
affected by residual positive magnification bias, produc- 
ing a systematic trend compared to blank fields. This 
trend is indeed expected under simple considerations, as 
shown in Sect. 18.21 and Fig. El A positive magnification 
bias is expected when the slope of number counts (with 
the approximation a = —d{log n)/d{log L)) is a < —1 
in the magnitude and redshift domains considered, i.e. 
Niensedi > L) = N { > L) X (scc e.g. Broadhurst 

et al. f995). This is indeed the case within our H band 
limited sample, according to the simple assumptions given 
in Sect. 18. 2I and the shape of the observed LF. We have di- 
rectly accounted for magnification biases using the lensing 
models. However, an additional magnification bias could 
remain in our sample due to a systematic trend: up to a 
given limiting magnitude, we tend to detect the sources 
with the largest magnification factors, instead of (or in 
addition to) the intrinsically brighter sources, as it hap- 
pens in blank fields. This systematic trend could slightly 
modify the slope of the LF derived for high-z sources. It 
is difficult to correct for without a complete mock simula- 
tion, assuming a shape for the LF of background sources, 
and then statistically correcting for this additional bias. 
Given the error bars obtained for the LF in Sect. 18.41 this 
residual magnification bias should be a second order cor- 
rection for high-z sources. However, the same trend could 
exist for very faint intermediate-redshift interlopers, such 
as the extremely-faint source A1835#35 (Richard et al. 
2003). This effect is presently uncorrected in our sample 
for obvious reasons, but it could be responsible for part of 
the discrepancy, because intermediate-redshift interlopers 
should mainly affect the bright end of the LF. 

The standard calibration used to convert the ^1590 lu- 
minosity into SFR, which assumes a standard Salpeter 
IMF under equilibrium conditions (i.e. constant SF over 
timescales > 10®~^ yr), is not necessarily appropriate for 
objects at such early epochs. 



Table 3. Summary of results obtained for the comoving 
Star Formation Rate density (in Mq yr~^ Mpc~^) by in- 
tegrating the Luminosity Function down to 0.3 Ll^^ in 
different redshift ranges, also presented in Fig. ^1 cor- 
recting for lensing, incompleteness effects and expected 
fraction of spurious sources. 



Assumption 



[6 — 10] First-category 

[6 — 10] First and second-category 

[8 — 10] First-category 

[8 — 10] First and second-category 



SFR density 
2.7 10"^ 

7.4 10"^ 

3.5 10"^ 
1.1 10"^ 
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Fig. 15. Evolution of the comoving Star Formation Rate 
density as a function of redshift. Different approaches are 
used to derive an order of magnitude estimate for this 
quantity, reported in Table Results from other surveys, 
uncorrected for extinction, are compared to our upper 
limits taken at face value. Data are compiled from the 
CFRS (Lilly et al. 1995, filled triangles), Connolly et al. 
1997 (filled pentagons), LBG work from Steidel et al. 1999 
(open squares), Fontana et al. 2003 (open circles), Iwata 
et al. 2003 (cross), Bouwens et al. 2003a (filled squares), 
GOODS (Giavalisco et al. 2004, open stars), different ACS 
estimates from Bouwens et al. 2003b (filled stars) and 
Bouwens et al. 2005b (filled circle). We also report the 
value derived by Bouwens et al. (2004b) and Bouwens et 
al. (2005a) in the Ultra-Deep Field ("UDF" labels). Our 
results (filled red diamonds) are presented for both clus- 
ters in the redshift ranges [6 — 10] and [8 — 10] : these val- 
ues are obtained by integrating the Luminosity Function 
fit down to 1/1500 = 0.3 Ll^^. Solid lines refer to the first- 
category candidates only (Am/f < 0.4), whereas dotted 
lines correspond to first and second-category sources. 



8.6. Spectroscopic follow-up 

We have started the spectroscopic follow-up of our sam- 
ple of high-z candidates with ISAAC/VLT. This survey is 
presently ongoing, and the final conclusions will be pre- 
sented in a further paper. Results on these observing runs 
have been (partly) published in Richard et al. (2003) and 
Pello et al. (2004a) , as well as a first preliminary summary 
in Pello et al. (2004b). 

To search for faint emission lines, we have systemati- 
cally explored the 0.9-1.40 /im domain {SZ and J bands of 
ISAAC), where Lya should be located for objects within 
the 7 < z < 10.5 redshift interval. We intended to detect 
emission lines with intensities ranging between 10^^*" and 
a few 10~^^ ergs s~^ cm~^ with a spectral resolu- 

tion for the sky lines of R — 3100 corresponding to the 
instrumental 1" slit width. The fraction of spectral band 
lost due to strong OH sky emission lines is of the order 
of 30%. Slit configurations were set to optimize the acqui- 
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sition of a maximum number of targets per night, with 
priority given to first-category dropouts. Secondary tar- 
gets were only observed when aligned with first priority 
targets. 

Up to now, our spectroscopic survey with fSAAC has 
targeted 2 priority candidates in AC114, and 7 in Abell 
1835 (4 "first priority" targets and 3 secondary ones). 
From this sample of 9 targets, 2/3 of the objects observed 
display emission lines. 5 sources have clear emission lines 
detected, and another one is still to be confirmed. The 
distribution between first and second priority targets for 
spectroscopy does not fully coincide with the present clas- 
sification in Tables IC.2I and IC.3I because it was based 
on an earlier version of the image reduction and analy- 
sis. For instance, z and SZ images were not available for 
Abell 1835 at this epoch. In summary, from 6 first pri- 
ority spectroscopic targets observed in the two clusters, 
we have clearly confirmed one candidate (A1835#8, orig- 
inally named A1835-1916, see Sect. I7.4ll . which is found 
to be a puzzHng source; two z > 7 candidates show emis- 
sion lines still to be reconfirmed; one candidate is found to 
be a low-z contamination, and two of them do not show 
emission lines. From the 2/3 secondary targets display- 
ing emission lines, one is a possible z > 1 source to be 
confirmed, whereas the other one is a faint low-z galaxy 
(A1835#35, z = 1.67). According to these preUminary 
numbers, the efficiency of our survey could range between 
~ 30 and 50%, with interesting low-z by-products such as 
A1835#35 (Richard et al. 2003). A large majority of our 
high-z candidates still need to be confirmed, either by a 
redetection of a faint emission line, or by the non-detection 
of other lines expected at low-z. 

9. Summary and conclusions 

We have obtained deep JHK near-IR photometry of two 
well-known lensing clusters, A1835 and AC114, plus z and 
SZ imaging of A1835 with FORS and ISAAC at the VLT. 
Our photometric depth reached SZ ~ 25.6, J ~ 24.4, 
H ~ 23.5 and Ks ^ 23.3 (Vega system), in addition to 
a minimum magnification factor of 1 magnitude over half 
of the ISAAC field of view. These images, combined with 
existing data in various optical bands including HST im- 
ages, have been used to select galaxy candidates at very 
high redshift {z ^ 6 — 10). The candidates have been se- 
lected with the dropout technique and two-color selection 
criteria appropriate for high-z galaxies. 

From our H band selected sample we have identified 18 
(8) "first and second-category" optical dropouts in A1835 
(AC114) up to ffvcga ^ 23.9 (AB - 25.3, uncorrected for 
lensing). Second category is defined here as objects de- 
tected in > 2 near-IR bands, the best-detected sources 
being defined as first priority. Among them, 8(5) exhibit 
homogeneous SEDs compatible with star-forming galaxies 
at z > 7, and 5(1) are more likely intermediate- redshift 
EROs. In both fields we have also identified a few ad- 
ditional dropouts detected only in the H band ("third 
category" objects), which satisfy our photometric selec- 



tion criteria. We have estimated the fraction of spurious 
sources expected in the different filter combinations, and 
corrected all the relevant derived quantities for this effect. 

Typically our candidates are magnified by a factor of 
1.5 (~ 0.44 mags) to 10 (2.5 mags), with average (me- 
dian) values of the order of 6.5 — 7.9 (2.3 — 3.5) for the 
two clusters. All high-z candidates turn out to be fainter 
than iJvcga ^ 23, with typical effective (i.e. lensing cor- 
rected) magnitudes ffvcga ~ 24 - 25 {Hab ^ 25.4 - 26.4) 
and fainter in some cases. Assuming standard SFR(UV) 
conversion factors, which may however be questionable for 
galaxies of such presumably young age, the SFR is found 
to be between few units and ^ 20 M0 yr~^. Their UV 
restframe spectrum, measured by the H — K color, seems 
to be very blue - a trend also reported for other high-z 
galaxy samples (e.g. Papovich et al. 2004, Bouwens et al. 
2004b). 

Taking into account the gravitational lensing effects, 
sample incompleteness and expected spurious detections, 
a first attempt was made to constrain the density of star- 
forming galaxies present at z ~ 6 — 10 based on lensing 
data. Integrated quantities, such as LFs and SFR den- 
sities, are affected by non-negligible uncertainties due to 
the large correction factors applied to this sample. The LF 
measured for LBGs at z ~ 3 seems to be slightly higher 
(by ^ 0.5 dex) but still consistent with the LF derived 
for our sample. The turnover observed by Bouwens et al. 
(2005) towards the bright end relative to the z ^ 3 LF is 
not observed in this sample. 

We have also estimated an upper limit for the cosmic 
SFR density from these data. Our values in the z ~ 6 — 10 
domain arc higher than the estimates obtained in the 
NICMOS UDF, even when the most conservative assump- 
tions and corrections are applied to the data. This differ- 
ence is related to the bright end of the UV LF for our 
candidates, i.e. L1500 ^ 0.3L*^3. This systematic trend 
with respect to blank fields could be due to to field-to-field 
variance, a positive magnification bias from intermediate- 
redshift EROs, and/or residual contamination by spurious 
sources. Given the error bars, residual magnification bias 
should be negligible for high-z sources, but a population 
of faint intermediate-redshift interlopers affected by posi- 
tive magnification bias cannot be excluded. At least one of 
such intermediate sources was spectroscopically identified 
in the field of A1835 (Richard et al. 2003). 

According to our simulations, in agreement with the 
first photometric results presented in this paper, the use 
of lensing clusters as gravitational telescopes seems to sig- 
nificantly improve the survey efficiency of z > 6 galaxies 
compared to blank fields. 

Given the uncertainties involved in the candidate- 
selection process, and the faint fluxes observed for our 
photometric candidates, the present results are to be con- 
sidered as a first attempt to constraint the population of 
6 ^ ^ 10 star-forming galaxies using lensing clusters. 
The present results and conclusions have to be confirmed 
and improved. Spectroscopic follow-ups are underway to 
determine the efficiency of our selection technique, and the 
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contamination level by intermediate-redshift interlopers. 
Additional deep photometry in various bands are being 
secured with HST, IRAC/Spitzer, and from the ground 
to improve the SEDs characterization of the high-z candi- 
dates. Increasing the number of lensing fields with ultra- 
deep near-IR photometry is essential to obtain more accu- 
rate constraints on the abundance and physical properties 
of z > 7 starburst galaxies. 
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Appendix A: Further improvements in the data 
reduction 

The following steps were introduced in addition to the standard 
scheme to improve the data reduction; 

— Object masks created by XDIMSUM were not well-suited 
for our fields centered on lensing clusters, because they did 
not correctly take the bright extended haloes in the clus- 
ter core into account. We improved the sky-subtraction in 
these regions by applying a simple threshold above the sky 
background in our images to create the object masks. This 
procedure greatly reduced the contamination close to the 
bright objects. 

— About 35 % of the images taken with the Hawaii-Rockwell 
array suffered from bias residuals, that appear more pro- 
nounced at the bottom and middle region of the detector. 
Before combining the individual sky-subtracted frames, we 
removed these residuals by subtracting from each line its 
median, with a rejection of pixels flagged in the object 
mask. 

— In the case of AC 114 cluster, about 20% of the individual 
H band frames presented strong low-frequency background 
variations, due to imperfect sky-subtraction and possible 
contamination due to the proximity of the moon. We cor- 
rected part of these residuals using a bidimensional fit of 
the large-scale background for each frame, after rejection 
of all pixels flagged in our object mask. 

— Before combining the frames into a final stack, we applied 
weight values accounting for slight variations in quality 
during observations, in such a way that the best-quality 
images will have the highest weight. Weighting was opti- 
mized in order to improve the detectability of faint point- 
like sources ; we computed individual weights using the 
following relation : 

weight oc [ZP x var x s^)~^ 

where we computed the individual zero-point ZP and see- 
ing s from the magnitudes and FWHM of the 5 bright- 
est unsaturated stars located in our field. The local sky 
pixel-to-pixel variance var was derived through background 
statistics inside a small region free of objects. 

— To check for the final accuracy of our absolute photomet- 
ric calibration, we compared the theoretical and observed 
colors of several cluster elliptical galaxies for which we had 
available spectroscopy, after having reduced the images as 
described below, and seeing-matched them to the worst 
value. The empirical SED template compiled by Coleman, 
Wu and Weedman (1980) was used to derive the expected 
colors for elliptical galaxies at the cluster redshifts. In the 
single case of the J band image in Abell 1835, we cor- 
rected a 0.1 mag offset in the zero-point. We also checked 
that the optical to near-IR colors of the brightest elliptical 
galaxies were consistent with the theoretical expectations. 
We found that our final absolute photometry is accurate to 
about 0.05 mag throughout the wavelength domain. 

Appendix B: Sample completeness and 
false-positive detections 

In addition to the purely photometric completeness effect in 
the detection of near-infrared sources, another incompleteness 



factor comes directly from our non-detection criterion in the 
optical bands. We statistically expect, by integrating a normal 
distribution function above this level, to measure a la flux for 
16 % of all dropouts in each optical band, which we did not 
include in our current selection technique. Because there are 5 
such non-detection filters for each cluster, this gives a statisti- 
cal completeness factor of Csampie = (1 — 0.16)^ ~ 0.42, which 
is unrelated to the observed NIR magnitude of the objects. 
We applied this additionnal correction to the final sample of 
optical dropouts. 

Our detection scheme was optimized to identify faint 
sources which are only detected on the near-IR bands, i.e. a 
subsample of the images. To evaluate the fraction of spuri- 
ous detections expected in our photometric catalogs, we con- 
structed a special noise image for each cluster and each near- 
infrared band, where all astronomical sources were removed 
while keeping the same noise properties: we subtracted by pairs 
sequential images obtained with similar seeing conditions, and 
then coadded them using the same procedure described above. 
The result is an image with the same noise properties com- 
pared to the final stacks, affected only by small residuals at the 
location of the brightest sources. We used SExtractor to de- 
tect objects inside these noise images, with the "double-image" 
mode and the same detection parameters as for astronomical 
sources. After masking the regions around the brightest objects 
and galaxy haloes to prevent any detection of source residu- 
als (as for the astronomical images), the number of objects 
detected in these noise images was compared with the num- 
ber of high-z candidates blindly selected as optical-dropouts 
within the same region on the astronomical images, for differ- 
ent ranges of magnitude and categories of optical dropouts as 
defined in Sect. 16.21 H-band magnitude bins have been defined 
to include a similar number of spurious sources in each bin. A 
source was conservatively included as positive (spurious) detec- 
tion in this table when Sextractor in the double-image mode 
was able to measure a magnitude for this source above the 
same detection level as defined for our candidates in Sect. 16.21 
Results are summarized in Table IBTI for each cluster. 

The fraction of false objects obtained with this procedure 
is overestimated compared to real catalogs, since no attempt 
was made to manually correct for obvious spurious detections, 
whereas all the dropout candidates included in the final cata- 
logs were visually inspected by at least two different persons. 
Their false-detection probability is therefore reduced compared 
to Table rO 

Appendix C: Additional tests on the reliability of 
optical dropouts 

In addition to the simulations presented in Sect. 15.31 on the 
completeness and spurious detections expected on the refer- 
ence H band image and on the other near-IR images, we have 
performed additional tests on the reliability of optical-dropouts 
identified on the near-IR images. 

C.l. Pseudo-x^ images 

Detection pseudo-x^ images were created from individual J, 
H and Ks band images in the following way ; each image was 
normalized by the noise la image, weighted by the square root 
of the corresponding exposure-time maps, and then all images 
were registered and averaged together. We applied the same 
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SExtractor detection parameters to these images, and com- 
pared the detection results. This procedure has some obvious 
limitations, because the final stacks are not independent from 
the original //-band detection images. For Abell 1835, ~ 89% 
of the first and second-category sources presented in Table 
are re-detected in the pseudo-x^ image, with positions less 
than 1 pixel off with respect to original centroids, except #13 



and #15. For AC114, 75% of the first and second-category 
sources are re-detected in the pseudo-x^ image (all sources 
except #2 and #5). By definition, third-category sources are 
clearly detected only in the H band image. However, sources 
#21 and #35 in Abell 1835 are also re-detected in the corre- 
sponding pseudo-x^ images. On the other hand, 60(86)% of the 
first and second-category high-z sources in Abell 1835 (AC 114) 
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Table B.l. Expected percentage of false-positive detections in our samples of third-category (only detected in H- 
band) and second/first-category (detected in H- and at least another near-IR band) dropouts for each cluster, as a 
function of the detection filters and the H-hand magnitude. The typical errors in these values are 15 and ^ 30 for 
Abell 1835 and AC114, respectively. 





3rd Cat. 






First / Second Category dropouts 






H range 


H 


SZ+H 


J+H 


H+K SZ+J+H J+H+K SZ+H+K 


SZ+J+H+K 


[mags] 






% of spurious detections in Abell 1835 (AC114) 






22.75 - 23.00 


0(0) 





0(0) 


0(0) 0(0) 








23.00 - 23.30 


35 ( 57) 


27 


12 ( 25) 


25 ( 33) 12 12 ( 33) 


12 


12 


23.30 - 23.75 


98 (100) 


66 


65 ( 88) 


61 ( 61) 28 27 ( 38) 


30 


12 


23.75 - 24.00 


100 (100) 


56 


96 (100) 


73 (100) 27 28 (100) 


40 


11 



(Tables EH and 01 are re-detected in at least another near- 
IR filter with a mean S/N> 2 within a 1.5"-diameter aperture. 
Thus, a majority of our optical dropouts, at least the first and 
second-category ones, have a significant signal detected in sev- 
eral near-IR bands. 

C.2. Transient objects 

We have estimated the contamination level expected in our 
images due to known transients, such as TNOs or supernovae, 
taking into account the combination schemes used to build 
the final stack, which typically reject 10 — 20% of the brightest 
pixels at this stage. The typical motions observed and expected 
for TNOs range between ~ (1 - 10)" /h (Trujillo et al. 2001). 
In the detection H band a TNO moves at least between 3 and 
6 pixels in 10% of the exposure time, depending on overhead 
details and the distribution of exposures during the observing 
period, thus making a TNO selection highly unlikely. On the 
contrary, supernovae events within the cluster, and particularly 
a type I supernova in the halo of a cluster member, cannot be 
removed with the usual scheme. We do not expect more than 
1-2 events per cluster in these deep observations, according to 
previous findings (e.g.. Gal- Yam et al. 2002). Very rare and 
unique events, such as lensed supernovae (e.g. Marri & Ferrara 
1998, Marri et al. 2000), or a tidal disruption of a star by a 
BH (cf. Stern et al. 2004) cannot be excluded either, but again 
such events are not expected to dominate the sample. 

C.3. Photometric stability 

The Abell 1835 data in the SZ band were obtained at two 
different epochs: 19 April 2004 (~ 4/i exposure, 77 images) 
and 15 May 2004 (~ 2h exposure, 45 images). Both series of 
images have identical seeing and photometric zeropoints. We 
used these data to check the compatibility of our magnitude 
measurements at different epochs, and particularly to set con- 
straints on the possible variability of faint optical dropouts in 
these bands. Among the 9 sources detected in the SZ band, 7 
of them show magnitudes and la error-bars (with asz ~ 0.1 
to 0.5 mag) fully consistent between the two epochs and with 
the final stack. A bright and clearly variable source (A1835#4) 
and a particular case (A1835#8), both discussed in Sect. EH 
are clear exceptions. 

C.4. Stacked images 

We have generated stacked J, H and Ks images for all first 
to third priority candidates, excluding EROs. A 10" x 10" re- 



gion has been selected around the ff-band centroid for each 
z>6 candidate in Tables 1021 and lOl for aU JHKs images. 
Multiplicative corrections have been applied to properly com- 
bine images coming from the two clusters with different photo- 
metric zero-points. After background subtraction, images have 
been averaged using IRAF routines and different pixel rejection 
schemes in order to obtain a "clean" region around the stacked 
source, although the final photometry does not strongly depend 
on the combination procedure. Representative results are dis- 
played in Fig. IB. II for first to third-category sources combined 
in different ways. We used Sextractor to measure the corre- 
sponding fluxes, best magnitudes and errors of these composite 
sources in the different bands and samples. Colors are obtained 
using the same apertures as for individual sources. Results are 
summarized in Table IUtI The S/N of the composite sources 
in all filters increases with the number of stacked images, as 
expected if a significant signal was present in a majority of sub- 
images and bands. Although these optical dropouts are likely 
to constitute a non- homogeneous sample, all the stacked series 
display a break between J and H (typically J — H ~ 1.5 to 2.50 
in the Vega system), and relatively blue H — Ks colors, thus 
a photometric SED corresponding to a dominant population 
of 2: > 8 sources (see Fig. or a noticeable contamination 
by spurious sources. The flux detected in the J band is clearly 
higher for the first and first -I- second-category dropouts, i.e. 
the brightest sources in the H band, as expected if a fraction 
of these sources are at 2 < 8. This result is the same when 
combining all dropouts at 2 < 8 (6 sources, all of them first 
or second-category), but the final S/N in this band is higher 
when blindly combining all first -I- second-category dropouts. 
The same comments stand for the Ks band, for which the best 
S/N is achieved for the first and second-category dropouts. 
The profiles obtained for the composite sources in the differ- 
ent fllters are all compatible with the seeing values in the H 
band (~ 0.5"), and slightly broader in J and Ks (~ 0.55" and 
0.45"respectively). 

False or transient sources stacked in these images will tend 
to enhance the trend towards artificially "blue" H — Ks and 
"red" J — H colors, thus providing an independent method to 
estimate the contamination levels reported in Table IHtI and 
applied throughout the paper. For 2 > 6 candidates stacked 
here, H — K colors provide a rough estimate of the restframe 
UV slope /3. For young starbursts, fi is found to be ~ —2.5 
in the local universe (Heckman et al. 1998), /3 ~ —3 for the 
2 ~ 7 — 8 candidates reported by Bouwens et al. (2004b), and 
it could reach values up to /3 ~ —3 for young starbursts at 
very low metallicity (see Fig. 1 from Schaerer & Pello 2005). 
As discussed in Sect. 18.11 our high-2 candidates tend to be 
extremely blue in H — K, but none of them could be excluded 
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from the sample to a 3a level on the basis of a clear /3 < 
—3.5, although two first category dropouts are close to this 
limit (A1835#7 and AC114#2, see also Sect.|0- Note that 
the contamination here refers to the integrated light in the H 
band, whereas Table Ib!T1 reports percentages in the number of 
sources. The results are the following: 

— The extremely blue colors of the third-category sample ar- 
gue for a significant contamination, between 60 and 65% 
of the total light for j3 ranging between —3.5 and —2.5 re- 
spectively. The faint population of third category dropouts 
could be dominated by spurious sources, as already ex- 
pected from Table IBTI 

— The second-category sample is compatible with a very low 
contamination level by spurious sources. It is negligible for 
/3 ~ —3.5, and it could reach up to 10% for (3 ~ —2.5. 

— The combined first-category image is clearly dominated by 
the brightest sources, and among them the two dropouts 
reach close to /3 < —3.5 at 2-3ct level (Sect.|H3l. This gives 
a contamination level ranging between 62 and 70% for the 
whole sample, whereas it is ^10% when these two sources 
are removed. Note that these sources are detected in several 
bands. 

— The combined first -f- second category (with all sources 
included) is qualitatively in good agreement with expecta- 
tions in Table lHrl a.e. 33 to 50% maximum contamination 
on the whole sample). 

Overall, the contamination levels based on UV slope con- 
siderations are optimistic in general as compared to the blind 
corrections applied according to Table IHtI 

C.5. Spatial distribution of optical dropouts 

We carried out a Kolmogorov-Smirnov (hereafter K-S) monodi- 
mensionnal test (Peacock, 1983) on the observed radial dis- 
tance of the sources with respect to the cluster central galaxy, 
because high redshift images are expected to be preferentially 
found around the critical lines for positive magnification bias 
(Broadhurst et al. 1995), and thus inconsistent with a homo- 
geneous distribution. Such positive bias is expected to be the 
case when probing the steep part of the luminosity function 
(see Sect. IHH and Fig. When comparing our candidates 
with a uniform distribution (after masking the usual regions of 
the image), we obtain an average to high resulting probability 
(~ 40% for AbeU 1835, ~ 2% for AC114). Thus, the K-S test 
indicates that the spatial distribution of our candidates has a 
low probability to be drawn from a homogeneous distribution. 
However, applying a similar K-S test for the spurious sources 
used in Sect. IKl^ gives us much lower probabilities (< 0.01%) 
for the false-positive detections to be drawn from a homoge- 
neous distribution. Such a high clustering level is expected for 
spurious detections, usually concentrated in specific regions of 
the image having a high noise level. When comparing the distri- 
bution of candidates and false-positive detections in the same 
K-S test, we find that these two samples are likely to have a 
different distribution (with probability < 0.01% to be drawn 
from the same parent distribution). 



J. Richard et al.: Abundance of « 6 — 10 galaxies from lensing clusters, Online Material p 6 

Table C.l. Photometric magnitudes and colors derived for the composite z>& source candidates in Abell 1835 and 
AC114 (see text for details). Upper limits correspond to 3 ct. 



Combination 


J 




H 




Ks 


J-H 


H -Ks 


1st category dropouts 


25.28 ± 


0.20 


23.68 


± 


0.09 


>24.7 


1.71 


< -1.0 


2nd category dropouts 


25.69 ± 


0.23 


24.02 


± 


0.08 


23.80 ± 0.13 


2.20 


0.48 


3rd category dropouts 


>25. 


6 


23.79 


± 


0.10 


>24.7 


> 1.82 


< -0.92 


1st + 2nd category dropouts 


25.38 ± 


0.17 


23.81 


± 


0.08 


24.22 ± 0.15 


1.53 


-0.32 


1st + 2nd + 3rd category dropouts 


26.18 ± 


0.31 


23.84 


± 


0.07 


24.76 ± 0.24 


2.45 


-0.86 
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Fig. C.l. (2 pages) Close-up of the best candidates in Abell 1835, showing the object and their surrounding 10 x 
10 arcsecs field. Objects satisfying the ERO criterion {R — Ks > 5.6) are presented above, and other first or well- 
detected second-category candidates are given in the next page. The FORS-Z band (non-detection criterion), and the 
ISAAC-NIR bands are displayed in linear scale, from —3a to 6a levels. Contours are for background level + 1, 2, 3, ... 
a respectively. On the left is the SED in the RzSZJHK bands and the best photometric redshift solutions obtained. 
When two possible solutions coexist, the higher redshift fit is displayed with a red dotted line. Fluxes values are given 
in fx, with units multiple of 10^^^ ergs s^^ cm^'^. 
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J. Richard et al.: Abundance of 2 ~ 6 — 10 galaxies from lensing clusters, Online Material p 9 




Fig. C.2. Same figure as lC.ll for best-detected first and second-category candidates in AC114. Close-ups correspond 
to the HST-i? and ISAAC-NIR bands. 



ID RA DEC SZ J H K <p z z range z /.liq fi -t'lsoo SFR notes 
(14:) (02:) 10'^"erg s"^ A/q yr'^ 



First-category dropouts 



#1 


0:58.278 


50 


26.65 


24.56 ± 0.18 


23.14 ± 0.11 


22.22 ± 0.10 


21.10 ± 0.04 


A 


[6.3-8.1] 


7.80 


1.27 


1.27 


1.28 


79.5 


83.5 


Ex. ERO 




















[1.3-1.9] 


1.74 


- 


- 


- 


- 


- 




#2 


0:57.538 


52 


49.85 


24.80 ± 0.22 


24.05 ± 0.27 


22.29 ± 0.11 


20.95 ± 0.03 


A 


[1.18-1.64] 


1.34 


- 


- 


- 


- 


- 


Ex. ERO (1) 


#3 


1:01.484 


51 


03.63 


24.03 ± 0.11 


24.54 ± 0.42 


22.69 ± 0.16 


21.71 ± 0.07 


A 


[6.7-7.3] 


7.05 


1.64 


1.66 


1.65 


12.7 


13.3 


Ex. ERO 




















[1.1-1.7] 


1.47 


- 


- 


- 


- 


- 




#4 


1:01.733 


51 


05.26 


24.31 ± 0.14 


23.50 ± 0.16 


22.82 ± 0.18 


21.90 ± 0.08 


A 


[6.2-7.8] 


7.53 


1.67 


1.7 


1.69 


39.6 


41.6 


Ex. ERO 




















[1.2-1.8] 


1.57 


- 


- 


- 


- 


- 




#5 


1:07.034 


51 


35.71 


25.82 ± 0.52 


> 25.60 


23.24 ± 0.28 


23.91 ± 0.55 


B 


[8.0-10.0] 


(9.00) 


1.92 


1.96 


1.96 


28.6 


30.1 




#7 


1:05.067 


50 


57.52 


25.81 ± 0.51 


25.34 ± 0.89 


23.39 ± 0.32 


> 24.70 


B 


[6.4-8.0] 


7.66 


1.53 


1.55 


1.54 


10.1 


10.6 




#8 


1:00.058 


52 


44.08 


25.36 ± 0.34 


24.99 ± 0.64 


23.40 ± 0.32 


24.00 ± 0.60 


B 


[6.1-8.0] 


7.38 


110.81 


86.53 


> 25 


< 0.6 


< 0.7 


Ex. (2) 


#10 


0:59.890 


50 


57.59 


24.18 ± 0.12 


23.74 ± 0.20 


23.45 ± 0.33 


21.77 ± 0.07 


A 


[6.8-7.7] 


7.33 


1.49 


1.51 


1.50 


33.1 


34.8 


Ex. ERO 




















[1.1-1.7] 


1.47 


- 


- 


- 


- 


- 




#11 


1:06.182 


50 


27.74 


> 26.90 


24.29 ± 0.33 


23.54 ± 0.36 


21.72 ± 0.07 


C 


[7.0-8.5] 


(7.75) 


1.31 


1.32 


1.32 


25.4 


26.7 


Ex. ERO 




















[0.8-2.5] 


1.78 














#13 


1:03.125 


51 


28.81 


> 26.90 


24.41 ± 0.38 


23.58 ± 0.38 


> 24.70 


A 


[8.0-9.2] 


8.29 


2.2 


2.26 


2.24 


14.9 


15.6 




#14 


1:04.209 


51 


54.55 


> 26.90 


24.77 ± 0.52 


23.63 ± 0.39 


24.53 ± 0.97 


B 


[7.3-10.0] 


8.54 


3.35 


3.5 


3.45 


9.7 


10.2 




#15 


1:02.540 


51 


12.84 


> 26.90 


25.40 ± 0.94 


23.63 ± 0.40 


> 24.70 


B 


[7.5-10.5] 


8.52 


1.81 


1.84 


1.84 


17.9 


18.8 




#16 


1:03.657 


52 


54.83 


> 26.90 


> 25.60 


23.64 ± 0.40 


23.08 ± 0.25 


C 


[8.0-10.0] 


(9.00) 


119.99 


55.04 


> 25 


< 1.6 


< 1.6 




#17 


1:05.013 


50 


27.11 


> 26.90 


> 25.60 


23.71 ± 0.43 


22.06 ± 0.10 


C 


[1.48-5.05] 


2.53 


- 


- 


- 


- 


- 


Ex. ERO 


#22 


1:02.551 


51 


30.06 


25.00 ± 0.24 


23.99 ± 0.25 


23.81 ± 0.47 


22.59 ± 0.16 


A 


[5.5-8.1] 


7.68 


2.28 


2.34 


2.31 


17.8 


18.7 




#23 


1:05.699 


51 


52.92 


> 26.90 


24.93 ± 0.61 


23.85 ± 0.48 


24.03 ± 0.61 


C 


[8.0-10.0] 


(9.00) 


2.61 


2.69 


2.67 


12.0 


12.7 




#24 


0:58.036 


51 


29.09 


> 26.90 


25.16 ± 0.75 


23.88 ± 0.50 


> 24.70 


C 


[8.0-10.0] 


(9.00) 


1.6 


1.62 


1.62 


19.3 


20.2 




#27 


1:04.299 


51 


57.19 


> 26.90 


> 25.60 


23.93 ± 0.53 


24.57 ± 1.01 


C 


[8.0-10.0] 


(9.00) 


3.54 


3.71 


3.68 


8.1 


8.5 




Secon 


d-catcgory dropouts 




























#6 


0:59.659 


50:54.73 


> 26.90 


> 25.60 


23.37 ± 0.31 


> 24.70 


B 


[8.0-10.0] 


(9.00) 


1.46 


1.47 


1.47 


33.8 


35.5 




#18 


0:58.890 


51:02.47 


> 26.90 


> 25.60 


23.72 ± 0.43 


> 24.70 


C 


[8.0-10.0] 


(9.00) 


1.47 


1.48 


1.49 


24.4 


25.6 




#19 


1:00.138 


52:05.20 


> 26.90 


> 25.60 


23.72 ± 0.43 


> 24.70 


C 


[8.0-10.0] 


(9.00) 


5.48 


6.03 


5.92 


6.1 


6.4 




#20 


0:58.860 


51:23.85 


> 26.90 


> 25.60 


23.72 ± 0.43 


> 24.70 


c 


[8.0-10.0] 


(9.00) 


1.67 


1.69 


1.69 


21.3 


22.4 




#21 


0:58.732 


51:53.86 


> 26.90 


> 25.60 


23.76 ± 0.44 


> 24.70 


c 


[8.0-10.0] 


(9.00) 


2.13 


2.18 


2.17 


16.1 


16.9 




#35 


1:00.693 


52:09.58 


> 26.90 


> 25.60 


24.00 ± 0.56 


24.25 ± 0.75 


c 




1.68 












(3) 



Table C.2. Photometric properties of Abell 1835 optical-dropouts. From left to right : identification number, astrometric position, near-IR photometry, 
photometric redshift quality (see Sect. 17^ for details), redshift range [zl — z2] and applied redshift z, magnifications at z = 6 (fie), at z = 10 (/iio), applied 
magnification (/i), unlensed luminosity at rest-frame A — 1500 A . 

Objects noted "ERO" satisfy the {R — K > 5.6) criterion for Extremely Red Objects. For EROs with a best-fit solution at high-z, the low- 2: solution is also provided. "Ex" is 
mentioned in the notes for objects excluded from the discussion, and underlined photometric entries correspond to forced undetections in a given near-IR band after visual 
inspection (see text for details). All undetections are displayed with photometric upper limits, using the 1 a limiting magnitudes inside a 1.5 " diameter aperture given in 
table [H 

Bracketed values of z are rough values based on the color-color selection diagrams. 

^ Near-IR counterpart of the sub-mm source SMIM J 14009+0252 (Ivison et al. 2000, Small et al. 2002). 
^ Also named A1835#1916 (Pello et al. 2004a) 
^ Also named A1835#2582 (Richard et al. 2003) 



N RA DEC 


J 


H 


K 


z 


z range 


z 


Me 






ilSOO 




SFR 


notes 


(22:) (-34:) 




















10'"' erg s^^ I 




Mq yr-^ 




First-category dropouts 


#1 58:49.777 46:54.95 


22.19 ± 0.05 


20.52 ± 0.02 


19.23 ± 0.00 


A 


[1.58-1.89] 


1.62 














Ex., ERO 


#2 58:49.040 47:21.94 


24.14 ± 0.32 


23.01 ± 0.22 


> 24.30 


B 


[6.3-9.2] 


6.90 


2.75 


2.86 


2.78 


10.2 




10.7 




#4 58:46.829 47:43.83 


> 25.50 


23.33 ± 0.30 


24.23 ± 0.99 


B 


[8.7-11.8] 


9.82 


13.68 


10.98 


11.20 


6.1 




6.4 




#5 58:46.505 47:25.96 


24.90 ± 0.64 


23.41 ± 0.32 


> 24.30 


B 


[6.0-10.0] 


8.51 


10.1 


12.62 


12.10 


3.3 




3.5 




Second-category dropouts 


#13 58:44.972 49:17.27 


> 25.50 


23.72 ± 0.43 


23.67 ± 0.59 


C 


[8.0-10.0] 


(9.00) 


1.54 


1.57 


1.57 


23.1 




24.2 




#14 58:53.511 48:37.85 


24.55 ± 0.46 


23.73 ± 0.44 


> 24.30 


A 


[5.2-9.6] 


6.91 


2.73 


2.84 


2.76 


7.1 




7.4 




#16 58:50.243 48:35.75 


24.55 ± 0.46 


23.90 ± 0.51 


> 24.30 


A 


[5.2-9.6] 


6.88 


23.21 


15.42 


21.35 


0.9 




0.9 




#19 58:43.844 47:35.31 


> 25.50 


23.98 ± 0.55 


23.69 ± 0.60 


C 


[8.0-10.0] 


(9.00) 


3.31 


3.51 


3.49 


8.1 




8.6 




Third-category dropouts 


#15 58:53.529 49:13.45 


> 25.50 


23.75 ± 0.45 


> 24.30 


C 


[8.0-10.0] 


(9.00) 


4.88 


5.31 


5.24 


6.7 




7.1 




#17 58:47.551 48:53.51 


> 25.50 


23.93 ± 0.53 


> 24.30 


C 


[8.0-10.0] 


(9.00) 


2.85 


2.97 


2.96 


10.1 




10.6 





Table C.3. Photometric properties of AC114 candidates. Caption is the same as in Table E2I 



